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This  report  was  prepared  by  Thomas  N.  Berne  is  ia  and  Robert  M*  Engle,  Jr.,  of  the 
Theoretical  Mechanics  Branch,  Structures  Division,  Air  Force  Flight  Dynamics  laboratory. 
The  work  was  conducted  in  house  under  Project  No.  1467,  “Structural  Analysis  Methods,” 
Task  Mo.  146702,  “Thermos last ie  Structural  Analysis  Methods,”  and  was  administered  by 
the  Air  Force  Slight  Dynamics  laboratory.  Research  and  Technology  Division,  Air  Force 
Systems  Command,  Wright- Patterson  Air  Force  Base,  Ohio.  Mr,  Robert  M,  Bader  is  the 
Project  Engineer  administering  Project  No.  1467, 


This  report  covers  research  conducted  from  July  19S4  to  July  1966,  The  manuscript  was 
released  by  the  authors  to  September  1866  for  pubJteatkm  as  a  technical  report. 

This  technical  report  baa  been  reviewed  and  is  approved. 
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ABSTRACT 


Tm  exact  solution  to  Zhs  problem  of  conduction  of  teat  tnroush  a  slab  is  develop 

formniaied  In  terms  of  m  infinite  series,  aUoraS^ SJi !£ 
mitf  boundary  conditions.  The  solution  is  programmed  is  FORTRAN  IV  for  er> 
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SECTION  I 

nnncr"CTiON 


T feo  conduction  of  hear  through  a  slab  is  governs*!  by  tte  following  partial  differential 
equation 

^  f  £T?  c>T 

37 1K  57j '  i’Cf~  !U 

For  cans  tent  thermal  dlifusivify,  this  eqtstfaa  simplifies  to 

*  7T  41-  !2) 

<JX*  S  df  : 

The  exact  solution  to tb:o «sfiat£oa  csn  be  formolafcd  iaierms  of  as  infinite  series.  This  report 
develops  its  asset  solution  for  arbitrary  initial  conditions  and  time  GepesSsnt  boundary  con¬ 
ditions.  The  solutioa  has  bees  programmed  is  FORTRAN  for  an  IBM  7S34  compter  srd  tbs 
scarce  program  listing  Is  contained  in  Appendix  t 


SECTIGN  H 

MATHEMATICAL  FGRMULATit>2i 


A.  BDUN&aBY  CONDITIONS 

like  general  soMton  of  Equation  (2)  mnst  satisfy  arbitrary  fn.tial  sM  time  dependant 
boundary  conditions  which  can  fee  expressed  in. ffe?  foltowjrg&rzn;  _ 


T  (xft)  =  f(x) 


t  =  <5 


(5! 


Ku:fr  +  K«- Ti  F<=4  (>-) 


x  =  o 


C4} 


K„-|^  *KaT=Ft  d,^  (X) 


x  =  L 


IS! 


1 


Is?  are  censteiite.Ssteotisg.dlSereatmiueg  of  these  coefficients  dictates  the  mode 

ov'hsaf  trsssfer  present  at  the  boundary.  By  various  combination®  of  constants,  she  imposition 
of  surface  temperature,  convection,  beat  fits  or  insulation  is  possible.  A  more  detailed 
diseus&fcss  oa  ths  interpretation  of  boundary  conditions  is  contained  in  Section  XU, 


-  dT 

if  nix****  k£ 

%  da  R*8 


*»■§** *  v 


t  in  i  &■ 


Figure  1.  Geometrical  Representation  , 

"The  &rml&iil&a  of  Duharaal’s  superposition  theorem  id  account  for  the  time  dependent 
'eonsfepr  conditions  necessitated  breaking  thl  sglutios.  into  a  transient  portion  satisfying 

a  steady  state  ptes  transient  wlfii  ^ro  inita&icondMons^ 

Since  tfe  restrtetfon  has  been  placed  on  the  bouaaary  conditions  Shut  the  time  dependency 
©as  be  exprasssd  uF  a  product  of  a  time  funefjosi  and  one  of  Us®  staadai-d  bouniary  conditions* 
"it  is  thprefo^e  possl^s  to  solve  the  eqnafi«ms  r%|ectfe^  the  time  variation  ^tad  then  modify 

Ijto  •;  .  ~'£Zm  'X-  '  -  '• 

-'iU  pf^bleSi  is  f|srai sitaplified  by  fe£®aMag  the  solution  into  two  parte  u  steady e  state 
{ s ,sat^fyin|  th8  arbitrary  boundary  cPadlilps, aad  atraji3iectportion,  ’I^.tJt,t  ) , 

satisfying  the  initial  temperature  distribution  apd  ifemogeneous  boundary  conditions. 

K:s?E:mzMmB  soisTm  l_  \r  '•  ;  .-  .,...•  A.  \  • 

■"  '.■  '  '  AT  -v:  '  '  - '  .  *-  -  -  . 

For  steady  ooadlfionsvvs  note  that^sOand  Ration  $$  simplifies  to 

\  ’  A  - /;  :  •-&  »  o  ; '  ■  ■•  -:  (s) 


The  polirtioa  of  itzza  &$ m&mte  Found  directly  by  Spte^c&ttoa  with  the  result; 

A/A-V  -  .  '  '  Ts  *  Ak*B 

Ws  now  bnpose  ifee  &?bifcm!y  iarmoary  conditions 


Ht  ■? 

S*  -- 


+  *«*  T  «  F0 


*  -  o 


O 
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£l 

d* 


Kg|  w  Kgg  T  -  Fj_ 


K  S  L 


m 


Substituting  Equation  (?)  into  (8)  and  (8)  yields 


K„  A  4-  K(iS  =  F9 


if  O’? 


K2,  A  +  Ka5liAL  +  13)  s  F, 


Cf  M 


The  Constanta  of  integration,  A  and  B,  can  now  be  evaluated  from  Equations  (10)  and  (it).  In 
order  to  keep  these  repressions  in  general  form,  the  solution  is  'Accomplished  by  Cramer's 
rule  with  the  result 


A  = 


B  r 


K12  fl  “ 


K  K22L  ~  (  K  »  K*« 
{Kgj’F  K2aL?  Fp  -  K,, 

^i2  ^22  ^“  ”  tKjjKj^^Kig  K2|) 


C*  TRANSIENT  SOLUTION 


(52) 


A  product  form  of  solution  is  assumed  for  liquation  (2).  and  designated  by  Xt  s),^(  t ),  ‘kib- 
stitution  into  Equation  (2)  yields  _  ' 


Rearranging 


XJi  $  =  ~X®' 


xtf 

X 


& 

c# 


CI4) 


(55) 


which  requires  that  each  of  these  functions  be  equivalent  to  some,  as  yetrarbitrary  constant. 

Then  setting  this  constant  equal  to  -/?  results  in  two  ordinary  differential  equations  of  tn« 
form 


+  <&  =  0 


(*S) 


X**  fi*  X-  0 


07} 


Equation  (8)  has  the  exponential- form  of  solution 


<P  =  C0  e  “S  ’ 

whereas  Equation  (17)  is  satisfied  by 

X(%)*  C,css  +  Cg  sin  /8x 


c»8j 


v!3) 


(20) 


The  solution  to  %f%ssA&w{$}  is  then 


?  (Si-U  av|o«  coe  fix.  *  Ca  sin  )9xj  [c0e“a^  *]  (20) 

This  immUitit  solution  mukiaxiisfy  the  initial  temperature  distribution  and  homogeneous 
boundary  conditions  m  mikxwm 


,  T  I 

«■  t)  -  f  (x) 

1  >0 

(3) 

| 

„  St 

i 

H 

* 

K.2^0 

X  =  0 

(2  i) 

•Os 

'i 

i 

i%4L  * 

r?X 

K22T=0 

x  =  L 

(22) 

J 

4 

"3 

-1 

-cd 

Bote  first  that  fee  constant  €?0  atm  be  eliminated  since  its  effect  can  be  included  in  and  C^. 
To  equate  the  remaps©  constants  substitute  Equation  (20)  into  Equations  (3),  (21)  and  (22). 

Substitution  of  Equation  (2G)  into  (21)  yields 


from  which  v/e  obtain 


«j8C3  +  K{,  $ 


C  -  K<»^  fs 

c*  “  T5T-C* 


(21) 


(2<3) 


At  this  point  it  becomes  necessary  to  impore  the  artificial  restriction  that  ^  0,  in  order 
that  calculations  performed  on  the  computer  remain  bounded. 

Substitution  of  Equation  (20)  into  (32)  yields 

[fCsfiC08 0L  -  Xa  sin  £t]c,  "*■  [ksi  0  gos.£l-5*  X22sfn  £l]c<i  ®  0  (25) 

To  obtain  a  nontrivial  solution  for  C,  and  CL,  the  determinant  of  their  coefficients  must  be  set 

12 

equal  to  aero.  This  yields  the  following  transcendental  equation, 

OL  z 


tan  2  * 


K*i*«r*“+  *2s  *«  *•* 


where 


and 


2  « 


^1)  ^2f~  ^lg 


(26) 


( 27 ) 


(28) 


--1 


1 
,  £ 

-  -i 
i 


1 


* 

._v 


'.4 
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Equation  <26)  has  infinitely  many  solutions  (eigenvalues),  and  we  shall  denote  these  by  ln  , 

where  n  -  0,l,E-**.  The  remaining  constant  Cg  is  evaluated  by  substituting  our  solution 

into  Equation  (3)  in  otder  that  the  initial  temperature  distribution  be  satisfied.  It  is  obvious  at 
this  point  that,  in  general,  arbitrary  functions  for  the  initial  temperature  distribution  can  rot 
be  satisfied  using  only  one  value  of  Z  and  C0,  We  are  thus  required  to  expana  our  initisi 

condition  and  our  solution  in  an  infinite  series.  Then  fee  coefficient  C2  becomes  CR  and  its 

evaluation  proceeds  as  follows,  The  total  solution  at  this  point  can  be  expressed 


T  U,(S  •  A*  +  B  +  S  Gn  [sin.  A,  »  cos  .So  *]* 


-«An£* 


tapering  the  problem  initial  condition  yields 


T  (x,o)  =  fix)  =  Ax+6  +  t0Cri  [sin  fth  x.  -  cci 


Kjf  £h 


o  -1 
a  j 


Bearrangiig 


f  (x)  **  (Ax+B)  -  jP  Cn  |isin\/3n 
n=o  L 


ft n 

x  -  — -77 — “*  cos 


ftn  *] 


Since  the  sices  and  cosines  form  a  complete  jset  of  orthogonal  functions,  the  C^a  can  be 
exmluated  by  midtiplyiigl»fe  sMss  o^E^atipa  ^1)  hy 

f  sin  ft}>  x  -  -  cos  IB#  x]  and  integrating  from  zero  to  L. 

Thus,  -  -  -  -  ....  -  . 


J  |f  (x)  ”  (Ax*i*B)  J  [sin  j3{ 


‘"Kin  " 


#*»*] 


/t-  ffl)  r 

X  cn  ftn  x  “ 

n=o  L 


K.I  fin 


ftn  4  fei  n  &M  x  -  cos  $N  xj  dx  (32) 


By  orthogonality  this  integration  produces  nontrivial  results  only  in  the  case  of  n  =  K. 
Therefore  /"  -  :  V  ::  : 


Ch  5 


/  [fx-(Ax4-e)l  [sin  ^  cos  ftn  xldx 

**0  .  -  ■g-  y  n[8  A 

[  (sln  ^n  «  ~  *“1^  co8  x]  dx 


5 
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The  denominator  of  CQ  c®  bs  evaluated  direct^  with  the  result* 


tb 


n  K!l^n 

„0  jl»  &  x““k7~c°5 


*  J—  J  f  f  K«^n  f  .  ,1 

d*s  IM  1“kT*)  +  *j 


-5-  sin  zn  cos  sr 


A  ']  !fl0 


.  2 

31  ft  ZR 


134) 


Collecting  the  formulations  required  for  problem  evaluation  leads  to  the  expression  of  the 
solution  for  time  independent  boundary  conditions  in  the  form 

'TUx,  t)  =  As  +  B 


where 


a,  L  [««>-(A.+B)}  [vn(«l]-d«  . afi% t 


■/»  tT' 


h  i*>] 


rK1£  FL-K«F0- 


^  “1  j 

f  *  ^Er'^kgaP^)  ~ 

3St  Kj,K8tL'^  -  j 


Y,  M 


=  [sin  Aj 


«  - 


% 

XT 


C09 


SJ 


135} 


mt 


-itn 


lot,) 


and 


fan  -  ZR  =  >] 


M?  ~  Z»/jL 

:  HZ, 


fi  -  K2IK„  Z|  +  K«.K,afc* 


52  -«<2 


-  f?3y 


-  V40): 


D.  MODIFICATION  OF  SOLUTION  FOR  TIME  DEPENDENT  BOUNDARY  €ONDIT*C‘SS 

Duhamers  superposition  integral  is  new  applied  to  the  solution.  Equation  (SS^rsO  a 
for  the  time  varying  boundary  conditions,  P.  3,  Hildebrand  give^taesGlutioa^ths^x^n* 

?{x,t}  ^  Ts{s,to)  f{P)i  (o)  4-  j  ^  m  d>,}  ^  fs> 

This  formulation  is  based  on  certain  limitations,  tov^rer,  which  must  be  eiinsis^fed.  T&e  first 
is  tbs  assumption  of  aero  initial  conditions.  This  restriction  is  eliminated  £ f  co^iaast^  a 

^HUSbian^  F,  B„  Introduction  fe>  Numerical  Analysis,  McGraw-ftSI  Booh  vorsipmy,  iso,. 
New  York,  1955. 
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xr^r&Kg  5  sto r 

to  bolt  iS^Sio^r^  *W3  result8  ”»  «“>  **<««.  Note  that  i(x)  =  o 

E.  SOLUTION  OF  INITIAL  CONDITION  PROBLEM 

BqStio!?  lifij'IS  of  ttetomogehM-astxjmda^K.omtions. 

A  -  0  t*  -  P  ¥zil“V  aA  sterfdy  elution  is  obtained  from  Equation  (?)  i  e 
result'  “  *  '  PI  -r  S  &"*  °°Mition  te  Equation  <35}  then  yfeL  the  dLtoed 


TfC  «M>  *L 


f  f{x}fyn  ix>]  dx 


fYn{st>] 


-a  ^ftt 


F.  UNSTEADY  STATE  SOLUTION 

lie  steady  state  solution,  Ts(x.® ),  employed  in  Equation  (35)  is  modifled  to  ’Cu.m)  A(t) 

• «s  as 

s  ,®)£(f)can  be  obtained  immediately  from  Equation  (7).  e-tog  boundary  coSton 
Equations  (4}  a sd  (5)  to  place  cf  (8)  azd  (9).  Tfe?  result  to 


where 


Ts  ^(U  ~  "2^  -  Glad'S 
.%Pl^LU?-KaFGA6(t) 


g  =  {K3*^agL?  Fo  &-W  "  PL  f>L  Ct) 

^12  K*j£L“D 

G*  TRANSIEKT  SOLUTION  (TIME  VARIABLE  BOUNDARY  CONDITIONS) 


„  ii . 

*«  dx  ■*■  K1£  T  s  F0  #3  (M 


K2f  3*  +  T  =  D 


122) 


AFFDL-TR-6&~109 


H,  COMPLETE  SOLUTION  FOR  THE  GENERAL  PROBLEM 

The  gsuaral  form  of  to  complete  soMlosa  «as  expressed  by  Equation  (41),  Collecting  the 
solutions  obtained  la  Equatiose  (42),  (43),  (49-  aad  (52),  and  substifestisg  into  Equation  (41), 
yields  to  ftol  result. 


f§P  M  = 


c  pi  s^n  X  % 

+  (  1 e)  +  J0  6  K  (M  dAj  *eft^  (s,t  ) 

r  rf  > 

+  {  Co  >  +  jfo  ®  ^SX)dX| f»T  ifqi)  4*  f  }C  (a,  * ) 


(53) 
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SECTION  W 


COMPUTER  PROGRAM  FOR  SERIES  TRANSIENT  . 
ANALYSIS  OF  SLAB  HEAT  TRANSFER  (STASH) 


A.  DESCRIPTION 

The  program  described  below  was  written  to  solve  for  the  temperature  distribution  in  a 
one-dimensional  rod  with  arbitrary  initial  conditions  and  time-varying  boundary  conditions. 
STASH  is  coded  in  FORTRAN  IV  for  the  IBM  7044-7094  H  Direct  Coupled-System.  Fifteen 
subprograms  malm  up  the  program,  each  of  which  has  a  specific  task  to  perform.  These 
subprograms  are  listed  below. 


MAIN 


reads  in  data,  sets  up  calculations^  and 
prints  the  results. 


SOLVE  i  -  solves  the  eigenvalue  equation  for  positive 
-  values  ox  DETK  (see  Appendix  II). 

SOLVE  2  -  solves  the  e^nvalue  eqyationfor  negative 

values  of  DETK.  ; 

SOLVE  3  -  solves  the  eigenvalue  equation  for  a  zero 

value  of  DETK. 

SC/LVE  4  -  solves  the  eigenvalue  equation  for  DETK  infinite. 
FHST  -  Simpson’s  rule  integration  routine. 

FUNCX.-  -  sets  up  the  integrand  for  the  x  integral 
FUNCT  -  sets  -up  tfieyiategrand  for  the  X.iniegr^f 
TABIN  -  reads  to  tabular  data,  if  present 
1NTERP  -  performs  linear  intexpolaii&s  ca  tabular  data 
PHIO 


PHD- 


dbfins3  ths  i’mm  varying  boundary  condition  at 

-  r 

-  (tefines  tto  time  varyii^  boundary  condition 
at  x  =  L 


PHXPRO  -  defines  the  derivative  of  the  time-varying 
boundary  condition  at  x  -  0  - 

PHJPHL  -  definasihs  derivative  of  the  time-?  varying 
boundary  condition  atx  -  L; 


FOFX 


-  defines  the  initial  conditions  in  the  rod. 


Figure  2  i3  a  simplified  flow  chart  depleting  tte  transfer  cd  Mdrmation  between  the 
subprograms  discussed  above.  .  * 


% 
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B.  INPUT 

Thwrs  are  basically  two  types  of  input  data  to  STASH,  the  physical  parameters  and  the 
problem  parameters.  Hie  physical  parameters  ate  the  characteristic*  of  tbs  rod  and  the 
conditions  to  which  it  is  subjected.  The  problem  parameters  are  the  accuracy  parameters,  Use 
calculation'  controls  and  the  print  controls. 

The  row  is  characterized  by  four  basic  quantises.  length,  mass  density,  thermal  conductivity 
and  specific  heat.  The  boundary  conditions  arc  Identified  by  the  indicators  K^,,  K^. 

The  magnitude  of  the  boundary  conditions  is  characterized  by  the  indicators  F^  and  F^.  Tbs 

magnitude  of  the  initial  conditions  is  similarly  characterized  by  the  indicator  Fx, 

Since  the  solution  may  be  required  at  any  station  on  2he  rod  for  any  given  pomi  in  time,  it  is 
convenient  to  specify  a  length  increment  and  a  time  increment  as  input  parameters.  A  final 
time  is  also  specified  to  terminate  calculations.  The  accuracy  of  the  solution  is  basically 
governed  by  three  factors,  the  accuracy  of  the  eigenvalues,  the  number  of  terms  in  tee  series 
portion  of  the  solution  and  the  aumner  of  intervals  takec  ia  fbs  numerical  integration  routine. 
In  the  interest  of  maximum  flexibility,  each  of  these  quantities  was  made  an  input  parameter. 

Information  may  be  input  to  the  programintwc  basic  forms.  The  first  t,pe  is  tin  date;  card. 
There  will  always  be  Severn  cards  inthe.data  deck.  If  the  tabular  data,  option  Is  used  there  may 
be  many  more.  The  second  type  of  input  consists  of  FORTRAN  IV  statements  ^duchmay  Jbc 
inserted  into  the  subprograms  defining  the  Initial  and  boundary  conditions  m  the  rod.  An 
example  problem  using  both  types  of  input  is  given  in  Section  V  C.  Detailed  instructions  ca 
inputting  the  data  cards  are  given  belowv  ln  the  following  format: 

(1)  Card  number  and  contents  : 

(2)  Program  name  for  contents 

(3)  :  Format  of  card  input  referenced  to  format  stetementmimber 

(4)  Description  of  each  variable  on  the  qard 

1.  Datacards,  -  '  :- 

Card  1  Intermediate  print  cations 

JDRINTfl),  JPRINT{2>,  JPflINT<3>  .  ..  "" 

5<m  FORiiAT{3D)  i 

JEWNTfl)  -  prints  series  portion  of  solution  term  by  term 

if  a  one  is  entered.  If  no  print  to  desired  enter  a  retro, 

JPjUN?f2)  -  prints  unsteady  state  portion  of  the  solution  if  e 
ems  to  festered.  If  no  print  Is  desired  ©iter  a  zero. 

JPRINT(S)  prints  solution  for  eigenvalues  if  a  one  is  entered, 

.  .  •  B’  imprint  is  desired  ester  a  sex©. 
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Card  2  Title  Card 

IUNIT,  TITLE 

1  FORMAT  (12,  I3AS) 

IUNIT  -  an  Indicator  which  prints  out  the  system  of  units 
to  fee  used  in  the  problem.  See  Table  I  for  a  list 
of  systems  presently  contained  in  the  program. 

TITLE  -  any  alphanumeric  Information  through  column  80. 

Card  3  Physical  Parameters  (all  must  be  in  consistent  units) 

L,  K,  RHO.  CP.  CELT  AX,  DELTAT,  T1MEF 

2  F0.1MAT  (7E10.0) 

L  length  of  the  rod 

K  thermal  conductivity 

RHG  mass  density 

CP  specific  heat 

DELTAX  increment  of  length  (100  increments  maximum) 

DELTAT  increment  of  time 

TIMEF  final  time  (Initial  time  is  zero) 

Card  4  Boundary  Condition  Indicators 
KX1»  K12,  K21P  K22 

3  FORMAT  (4E10.0) 


m 

indicator  for  at  x  =  Q 

K12 

indicator  for  T  at  x  =  0 

K21 

.IT 

indicator  for  at  x  =  L 

d* 

K22 

indicator  for  T  at  x  =  L 

Funotion.  Multiplying  Factors 

FO,  FL,  FX 

; 

4  FORMAT  (3E1G.0) 

FO 

coefficient  on  function  $0(t> 

FL 

coe.'  aient  on  function  $  ^(t) 

FX  coefficient  on  function  ffe) 


TABLE  I 


SYSTEMS  OF  UNITS  STORED  INTERNALLY 
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Card  6  Calculation  Parameters 

NTERMS,  KSTEPX,  NST£PTf  NTAS(l),  STAB  (2),  NTAB»3),  NTAB{4), 

NTAB('.i 

5  FORMAT  {315,  5X1) 

NTERMS  -  number  of  terms  in  fee  series  portion  cf  the 
solution  (100  maximum) 

NST3PX  -  number  of  intervals  for  The  x-iwtegr ation 
NS1EPT  -  number  of  intervals  for-  the  A  -integration 
NTAB{1)  .  -  flag  for  table  1  N 

NTAB(2)  -  flag  for  table  2.  NTAB(F}  =  0  Do  cot  use  Table 
NTAB{3)  -  flag  for  table  3  i 

%  -  NTAB(4)  -  flag  for  table  4  j  NTAB0)  =  1  Use  Table- 

NTAB<5)  -  flag  for  fable  5  } 

!f  NSTEPX  or  NSTEPT  ie  garo  fee  program  sets  fee  value  of  the  respective  integral  to  zero. 

Card  7  Eigenvalue  Solution  Parameters  : 

••  imk'pmmd.;  >'  . 

€T  FORMAT  ;M0.O*35) 

LIMIT  *•  difference  between  two  successive  iterations 
necessary  te  define  convergence  to  a  root, 

IT'EBMJC  -  ;  maximum  number  of  iterations  to  be  made  in 
seaychlag'fqr  sach  sigenv^us. 

If  no  tabular  data  is  to  bo  used,  this  isfee  last  caiid  ir,  the  cfete  deck.  If  tabular  data  is  to  b3 
an  input,  fat  wove r*.  she  following  format  twill  ba  used. 

Card  1  TaMe  Humber  and  Cvunraents 

NT  ABLE,  COMMENTS  ./ 

FORMAT  (I5„  20s,  &SH  )  ;  .. 

NT  ABLE  table  number 

COMMENTS  aay  alphsuumsrio  information  in  columns  ?jj  through  80. 
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Card  -2  Tabular  Data  {{2  to  50  data  cards  per  table) 

INDVAR,  DBPVAR,  COMMENTS 
FORMAT  (5X,  2E10.0,  55H  ) 

3SDY.4R  independent  variable 

DEPFAR  dependent  variable 

COMMENTS  any  alphanumeric  Information  in  columns  26  mrough  80 
Card  u  End  of  Table  ' 

>  N 

FORMAT  (15) 

N  \  negative  <ti  table  number 

There  are  fivd  tables  provided  n\  thcprogram,  which  are  assigned  as  follows: 
Table  1  PHIO  -\V  . 

Table  2  PtiiL  .  ■  • 

. /Table  3  PSIPRO  5  - 

'  .  Tablft:4  V  PHFTRL  i  '  -  ,/  •  ? 

Tables  FOFX 


Figure  3.  shows 
XV  C. 


d  symbolic  data  deck.  A  sample  problem  i3  generated  in  detail  in  Section 


Subprogram  input  Cards  \  ^  .  I  ?  .  1  - 

f  if  ths  tabular  data  optionis  notuse&STA^F"  2&.lcuUies  tbs  required  functions  internally  using 
FORTRAN  statements  as  loaded  in  the  subprograms  at  compilation  time.  The  affected  sub- 
programs  are:  .  -  .  - 

FUNCTION  RHIO 

FUNCTION  PHIL  ' 

FUNCTION  FH1FRO 
FUNCTION  PHIPRL 
FUNCTION  F$FX 

As  the  Initial  or  boundary  conditions  change,  cards  containing  the  iimGttoaal  statement  of  tbs 
variation  must-be  insertesa  Since  all  the  above  functions  have  a&  associated  multiplying  factor 


Figure  3.  Ii^iut  Data  Forrest 


13 
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it  is  convenient  to  use  normalised  functional  statements  in  the  subprogram,  Tims,  if  we 
entered  the  following  card  in  the  FUNCTION  FOFX 

r  FOFX  =  1.0 


we  would  obtain  the  initial  condition 


££*.}  =  FX 

The  program  (Appendix  I)  contains  subprogram  statements  corresponding  to  the  following 
initial  and  boundary  conditions. 

^>q{  i)  =  const 

<pL  (t)  =  const 

4>'0®  =  o 

$  (t)  -  o 

f(x)  =  const 

G.  SAMPLE  PROBLEM 

Consider  a  rod,  ten  inches  long,  having  the  following  properties: 

K  =  10  Btu/hr-ft-*F 

p  «  SOO  lb  /ft3 

r  XXL 

c  =  0,1  Btu/lb 

I»  m 

We  wish  to  obtain-a  time  history  of  the  tempera  ture  -  distribution  through  the  rod  which  Is 
subject  to  the  following  boundary  conditions 

T(G,  t)  =  t 

T(L,  tj=  0 

The  initial  conditions  imposed  on  the  rod  are: 

T(  X,  0)  a  0 

For  purposes  of  illustration  the  intermediate  print  option  will  be  called  out  on  card  ore. 

The  first  task  in  setting  up  the  problem  Is  to  decide  upon  a  system  of  units  to  employ.  Since 
the  length  of  the  rod  is  given  in  inches,  we  shall  choose  the  inch  as  the  unjt  of  length.  For  a 
transient  study  a  small  time  unit  is  desirable,  hence,  the  second  becomes  the  unit  of  time.  The 
remainder  of  ths  system  of  units  is  determinedly  the  temperature  and  mass  units.  Going  now 
to  Table  I,  we  obtain  the  correct  value  of  IIW1T.  This  indicator  along  with  a  suitable  title 
becomes  card  two  in  Figure  4. 


19 


APFDL-TR-66-103 


Card  three  contains  the  physical  parameters  of  the  system  in  tbs  system  of  units  called  lor 
by  tbs  indicator,  TUNIT.  Toss?  parameters  are: 

L  =  10  in.  K  -  0*002314  Btu/te.-sec-“F  RHO  »  0,2895  lbm/ito3 

CP  =  0.10  Biu/Ife  DELTAX  =  0.5  in.  DEI/TAT  -  100  sec  TIME/  =  1000  sec 
m 

Sines  we  have  only  tempeiafure  boundary  conditions  we  set  Kll  --  K21  -  0  and  K12  -  K22  -  1, 
This  information  Vs  entered  os  card  four. 

The  functions  defining  our  initial  conditions  are 

90(t)  -  t,  #L(t)  =  0,  ffc)  =  0,  ^{X)  =  ls  «^{X>  -  0 

We  can  take  advantage  of  the  functions  already  stored  In  the  p^ogr am  by  using  and  F _  to 
zero  out  the  appropriate  functions.  We  then  can  use  Fq  -  1  to  bring  in  the  other  boundary 
condition.  This  is  shown  on  card  fire. 

The  calculation  parameters  are  entered  on  card  six  A  series  solution  containing  twenty-five 
terms  is  completely  adequate  for  this  problem.  Since  the  initial  condition  function  is.  sesv,  we 
set  rSTEPX  equal  to  zero.  For  a  Simpson’s  rule  integration  eehes »  fifty  steps  shout*  suffix 
for  NSTEFT.  Our  choice  of  multiply  ing  factors  os  card  five  enabled  03  to  do  much  of  the 
function  calculation  Interssdly.  However,  to  eliminate  recompiling  any  portion  of  the  program 
we  used  tables  ibdaflBe  thefunetfons  tf>  (t)  and  <(>  (X),  Thus we  set  NTABflJ  and  HTAB{3) 
“  equal  to  Gse  and  the  rest  equal  to  zero. 

For  a  solution  with  temperature  boundaries  only,  the  eigenvalues  become  simply  mr  .  Thus, 
the  eigenvalue  parameters  have  little  meaning.  However*  for  a,  more  complex  case  they  would 
have  a  significant  effect  on  the  solution- so  these  pariimeters  should  be  made  as  stxingent  as 
required.  Typical  values  are  entered  on  jard  seven. 

Our  choice  not  to  recompile  any  subprograms  leads  to  the  use  of  tables  one  and  three.  The 
first  card  in  each  table  is  stable  designation  nsEiber,  The  following  cards  contain  data  points. 
The  last  card  of  each  table  contains  fee  G2gative  of  tbs  table  designation  r-umber  and  is  a  flag 
signalling  the  end  of  the  table. 

This,  then  is  the  data  deck  for  the  sample  problem.  The  assembled  deck  is  shown  in 
Figure  4. 

D.  RESTRICTIONS 


Certain  restrictions  mu3i  be  adhered  to  in  ordsr  for  tbs  solution  to  be  successful.  Violation 
of  these  restrictions  will  usually  produce  an  error  message  from  the  computer  program. 

a.  A  consistent  set  of  units  must  be  employed.  An  indicator  is  provided  on  the  title  card 
which  will  label  the  system  of  units  on  the  output.  If  this  indicator  is  omitted,  the  following 
error  massage  is  printed. 

SYSTEM  OF  UNITS  NOT  SPECIFIED.  1UNIT  NOT  ENTERED  OR  ZERO. 

This  message  merely  informs  Re  user  of  this  omission,  execution  ox  the  problem  is  not 
terminated. 
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b.  and  K^2  cannot  be  sere  simultaneously.  Similarly,  sad  cannotbe  sero 

simultaneously,  These  situations  lead  to  an  undefined  boundary.  The  error  message  below 
results  from  this  case. 

r  BOTH  INDICATORS  AT  ONE  BOUNDARY  ARE  ZERO 

c.  Kj  2  cannot  be  zero.  This  Is  a  somewhat  artificial  restriction  imposed  by  the  formulation 

of  the  problem,  fc  a  physical  sense  it  prevents  jq  possibility  of  the  uasolvable  Neuman 
problem.  If  is  zefo  ihs  following  error  message  is  printed. 

FORMULATION  BOX  NOT  PL  illT  KI2  TO  BE  ZERO. 

d.  The  number  o i  integration  intervals  muat  h>  iven.  This  restriction  arisesfrom  the 
computer  formulation  of  Simpson’s  rule.  If  an  ts&  v  .  r  umber  is  entered,  the  folfowicgerxor 


NUMBER  OF  INTEGRATION  INTERVALS  IS  NOT  EVEN 

e.  The  computer  program  generates  an  error  message  if  tce  time  increment  or  the  length 
increment  Is- zero  dr  iKgative.  This  message  reads:  :  . 

TIME  OR  LENGTH  INCREMENT  E  ZERO  OR  NEGATIVE 

X  The  initial  time  for  each  program  Id  zero,.  .  -  : 

:  g.  The  program  uses  even  increments  of  time  and  length.  The  maximum  number  of  length 
increments  ia  ona  hundred. 

E.  QUTJPUT  '  •  , 

The  output  generated  bySTASHconsists  of  two  segments;  the  input  dabt  display  and  tempera- 
tm>e  profiles  wMch  are  always  generated;  andidie  intermedia  print  vdiicb  is  controlled  bylbe 
fird'  eard  intfaedatadeck.  Aftor  readii^  vfee  date,  STASH  prints  it  out  along  with  suitable  titles 
and  headings  as  shown  in  Figure  b.  if  the  eigenvalue  solution  is  requested  a  table  of  the 
eigenvalues,  asditeratioiis,  is  printed  as  shewn  in  Figure  8.  Intermediate, print  options  giving 
the  valt&s  of  the  series  and  theunateady  state  produce  output  as  sbownin  Figure  Yior  each 
station  alorg  the  slab  at  each  time  step.  Figure  7  also  shows  the  fosan  ofthetemperature 
profiles  asg^narated  ateach  time  step.  . 


Flgxiro  Simple  Output  (C&ce  4)  Showfrg  Dlspliy  oHnpa, 
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SOLUTION  FOR  EIGENVALUES 


ROOT  NO 
1 
2 

3 

4 

5 

6 

7 

8 
9 

1$ 


0.157G7962E  01 
0.47123838E  01 
0.78539814E  01 
0.10835574E  02 
0* 1413710SE  02 
0.172787S9E  02 
0t  204203525  02 
0, 235S1S44E  02 
0.23703537E  02 
0.23845123E  02 


Figuxe  S.  Sample  Outpat^Qase  4}  Showing  Eigeavsilue  Solution 
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SECTION  V 
CONCLUSIONS 


Several  clasps  of  problems  were  run  to  cheek  out  fee  program.  The  results  were  compared 
with  a  finito“dAfference  heat  transfer  program  (LTA)  which  was  developed  by  Lockheed 
Aircraft  Corporation*  These  results  were  examined  for  accuracy  and  speed  of  convergence. 

The  intermediate  print  feature  of  the  program  was  used  to  determine  fee  convergence, 
Examination  showed  that  convergence  was  quite  rapid,  usually  less  than  ten  terms,  for  constant 
ocundary  conditions,  away  from  time  zero.  More  terms  were  needed  in  the  vicinity  of  zero 
time  to  produce  convergence.  For  time- varying  boundary  conditions,  the  series  does  not 
convene  very  quickly.  However,  a  study  of  the  solution  convergence  showed  that  a  twenty 
term  series  using  one  hundred  integration  steps  obtained  results  within  one  percent  of  the 
LTA  solution  for  times  exceeding  one  hundred  seconds.  At  smaller  times  &  fifty  term  series 
with  one  hundred  integration  steps  was  required. 

Appendix  HI  contains  fee  results  of  five  check  problems  compared  with  the  results  from  fee 
LTA  finite  difference  program. 

No  comparison  is  mad©  at  zero  time  since  fee  program  obtains  these  values  from  the  initial 
conditions  rather  than  from  a  series  calculation.  The  curves  are  plotted  from  data  generated  by 
the  program.  Gas©  I  is  the  sample  problem  detailed  in  Section  XV  C.  All  other  cases  used  fee 
same  physical  parameters. 
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APPENDIX  I 

COMPUTER  PROGRAM  SOURCE  LISTING 


non  0000000000 


AFFDL~TH~8S~1G8 


eiSFTC 


STASH 

MAIN 


PA  P 

M94/2#XR7 


1 

2 

3 

4. 

5 

tf 

15, 

2C 


STASHOOl 

SfASHC'02 

SfasKOci 

SrA5f«J04 

St^HOOS 

ST&SHaUe! 

STA5H0O7 


A  GENERAL  SGUJTICH  TO  THE  ONE-D 5 PENSIONS.!  Mt'41  T&W$£Eg  PSri&LE 
t*  I i’H  TIHE-OEPENOENT  80UNDARY  CONDITIONS  AND  AS8  USAS Y  £N  iflAL 
CONDITION 

VERSION  2  STA5H0O7 

VERSION  2  CALCULATES  ZERO  TIRE  USING  BOUNDARY  ANJJ  fN!  f«L  CfWiHT 10STASH005 
VERSION  2  INCORPORATES  A  MULTIPLIES  ON  THE  IN m*t  CONSr<TlON  HJNCTSTASHOGO 
AND  BOUNDARY  CCNOITION  FUNCTIONS  TO  PERSjiT  THE  USE  OF  NORMALIZED  STA'^-XUo 
FUNCTIONS  c  SlVlSHOll 

TFUNCpXFUNC 

TEMPI  1005,  TITLE!  1 3)  »EIGENUOS) 

JPKINTI3J 
NTABI5J 


EXTERNAL 
DIMENSION 
DIMENSION 

DIMENSION  _ 

REAL  tai8N«.KUBN2 
REAL  K11K2I >KI2K22 
SEAL 


KfcAS.  pKit'KH 

SEAL  L!HIT»t*K,lAM8DA»NlJHXS$eNUNSSeKFPHI,KFMHI0,*FPm\*Kil»>a2 
T.*K2l»K22t  N»KTERHl,KTERM2 


COMMON 


STASKQ12 
StASHOl* 
ST AS HOI 4 
STRSKQU 
^TASH0T6 
STASH&IT 
STASH018 
S7ASH019 
^-/ASH020 


IX,0cLTATpXlltK12>K21,K22pFX,NTAS  r 

SlMrkiiZ 

COMMON/ftOCTS/NTESPSjL* !TERHXpKilK2I «KI2K22* LIMIT 

5TASH023 

COMMON/PR INT/JPRINT 

STASH024 

STASH025 

FORMAT  STATEMENTS 

STA5H026 

STASH027 

FORMAT? I2» 13A6? 

STASH028 

FORMAT! 7E 10.0! 

STASH029 

FCRHAH4E 10. CJ 

STASHCUO 

FORMAT! 3E 10 ,0) 

STASH031 

rOft«ATC3I55in 

STASH032 

F0RHAT!E10.0»I51 

STASK013 

FORMAT! lHl*29Xp 13A6) 

STASH034 

22 

30 

31 


FGRHATIIH*p55X, I 9KNULTX PLYING  EACY0R5/6+X,  3KF0R/45X»  40H80UNDARY  ANSTASH04I 
15  INITIAL  CONOIT ION  FUNC7I0NS//30X* 4HF(0J»30X,4HFCU ,30X,4HFIX)//ZSTASH042 
24X*2!£t5.8*l9X}pE15.8l  STASH043 

FORMAT! IHQt29X»4(EI5o8»AX} J  STASH044 

FORMATI IHA,SlXf30HEI6ENVALUE  SOLUTION  PARAHET£RS//42X»8HACCURACYp  3STASH045 
•OX »9HNL‘M8ER  0F/44Xp5HLIKIT»31X, 10H£T£RATIOMS//39X,El!>.&,26X,!5}  STASK046 
FORMAT! 28X*4HINCH» 14X»4HSLUGp 14X, 3HSEC, 16X» SHPOUNOp 13X, 10UFAHRENHEST4SH047 
ITS  STASH04S 


13 
66 

120 

UTS  STASH04S 

130  F0RHAT!28X,4HINCHf14X,4HSLUG,14X,3HMIN*16X,5HP0UNDfI3Xf IOHFAHRENHESTASH049 
M-T-S  STASH05Q 

140  FORMAT !28X»4HiNCHf  I4X*4HSLUG»  I-5X*2HHR»-t6X»5HP0UND»13X»  10HF  AHRENHEI STA5H051 
l?  I  STASHD52 

150 


l&O 

17S 

ISO 

190 


IT  I  STASH052 

FORMAT!  28  X,4HFCUT*14X ->4HSLUG*  14  Xp  3HSEC »  16Xp SHPOUND*  1 3X»  lOrSFAHRENHcSf  A5H053 
IIT)  STASM054 

FQRMATr28X»4KFC0T*l4Xfr4HSLUG*l4Xp3HHlN»I6X»5HPOUN0t13X,IOHFAHR6NHESTAS?J055 
ITS  STASH056 


1!T1  STASH056 
FORMAT  USX, 4HPQQT*l4X«.4HStUGtl5X,2HHR*l6X,5HPOUN0,i3X,10HFAHRS«HE!STASHO5? 
IT)  STA5H068 
FORMA Tf 28 X p4H! NC Hp 14X» SNPOUNDp 1 3X#3HSEC» I6X, 5HPQUND, J3X, 10HFAHREMHSTA5H059 
1EITI  STASH060 
FORMAT!  2aX,4HA?-iCH»14Xp5HPOUN0,13X,)3HMIN,16X,5HPGUN0pl3Xpl0HFAHB6NHSTASKO&i 
IE JTI  STASH062 


28 


UOU  U  (->  O  O  O  U  O  O  O  U  O  <J  o  o  o 


A?FDL*-TR“6tt-108 


303  F0ftMAT(  JHS,<i*y,22HCAt.CULATI0^  PARAMETERS//26X v  i  5HNURBER  OF  TERMS*  1STASHG63 
!5X*HMNUMPt»  OF  INTERVALS.  I  JX*19MNUMBER  IF  iN!£RVALS/25X»i6!TlN  FHE5TA5N064 
2  SUMMATION, 14X.21FFQR  THE  X  INTEGRAT'Gn, 11X,?1HFGR  THF  T  INTEGRA! I  ST ASH065 
3CN//30X,  i5,.?BX,l5,26X»I5)  S1ASH066 

333  FGRH*.TI IHA.6TX  .SHA1PHAJ  STASH067 

1002  FCRHAT126X5S£).5.B)  STASH068 

1003  FORHAm.X*ElS,e}  STASH069 

1100  FORMAT { 23X,4HINCH,  14X,  5HPQL<N0»14X,  2HHR, !6Xt5HP0UND,  i  3X» 10HFAHRENHEST  ASHO  TO 

1IT)  STASH071 

IliO  FORMAT { 28X,4HrCGT« I  AX. 5HP0UND, 1 3X,  3HSEC, 16XV5HPQUND, 13X, IOHFAHRENHST ASHQ 72 
IE  IT)  STASH0  73 

1120  FORMA  n  28  X.4HFC0T, I4X, 5HP0UNG, 1 3X#3HHIN, 16X,5HP0UND» 1 3X» IQHF AHR6NHSTASH074 
IE  IT)  STASH075 

SI  30  FORMAT! 28  X*  4HFCQT t  14X, ShPOUNO,14X,2HHR, 16X,  5HPQUND* !3X* iQHPAHRENHESTASHQ  76 
HTJ  STASH077 

1U0  FORMAT! lHA,4SX.4GHPHySICAL  CONSTANTS  TO  DEFINE  THE  PROBLEK//59X,i5$TASH078 
1F-SYSTEM  OF  UNI TS//27X , 6HLENGTH* 1 3X»4HMAS5» 14X, 4HTINE , !4X*6HWE IGHT,S7ASH079 
213X,11HTEMPFRATURE)  STASH080 

2002  FORMAT{7Xf4hTIPE,22X,42HTcMPERATURE  DISTRIBUTION  ALONG  THE  ROD  AT  STASHOSl 
1  «r5. 3* ION  INTERVALS)  STASH082 

3003  FORMAT! 1H2)  STASH083 

3333  FORMAT! 1HQ* 58X*E 15.8)  STA5H084 

4 OCA  FORMAT (lHAf63X*4HCETK)  STASH085 

4444  FORMATS 1H0,58X,£ 15*8)  STASH086 

5000  FORMAT { 51 1 )  STASH087 

6100  FORMAT;iHl!>25X,2HTa*EI5.8,2UA>2HX*JEl5.8//45Xf26H?.i:RItS  PORTION  0FSTASH088 
S12LUTI0N//14X, IHN»I6X,4HTERX,19X*9HSUMNATI0N//I  STASH089 

6110  FORMAT {I0XtI5.10XjEl5.S»I0Xv£15«8)  STASH090 

62 CO  FORMAT! IHO, 40X, 34HUNSTE AOY  STATE  PORTION  OF  50LUTIGN/235X, IHX,25X»STASH091 
•4HTEHP/T)  STASH092 

6210  F0RMAT{28X,f 15.8 ,2QX»Ei5.8)  STASH093 

I ERROR-O  STASN094 

STASH095 

INTERMEDIATE  PRINT  OPTIONS  STASH096 

STASH097 

REAQI 5,5000)  JPRINT  STASH098 

STASH099 

JPRINT=i, PRINT  INTERMEDIATE  CALCULATIONS.  JPRINT=0,00  NOT  PRINT.  STASH100 
JPRINT »  l  }•*?  GRIP'S  PORTION  OF  THE  SOLUTION  TERM  BY  TERM  STAJHIOl 

JPRINT! 2) -UNSTEADY  STATE  PORTION  OF  THE  SOLUTION  STASH102 

JPRlNrnj-SOtUTiCN  for  THE  EIGENVALUES  STASH103 


READ  I5,1)IUNIT, l T I fLEII) , 13) 

R6ADI 5,2) L»K,RFO,CP»DELTAX»PELTATiT IMEF 
READ! 5,3)Kll,Xl2,K2l^K22 
REA0;S,4)F0,FL*FX 

READ! 5,5)  NTERMS «NSTEPX,NSTEPT *NTAB 
RfcAOI 5,6)  LINI T* ITfRMX 


PRINT  INPUT  DATA 
TITLE 

WRITE  ?6,itM  TITLE (I),1*U  13) 

SYSTEM  Or  UNITS 

WRITE  f  6* 1 140) 

IF{ IUNIT.EO.O)  GO  TO  5500 

GO  TO  1 102, 103, 1 24, 105, SG6, 107, 108, 109* XIO, 11 1, 112, 1 13} *fUNIT 
102  fcRIT6(6,120) 


STASHI01 
STASH102 
STASH103 
ST4SH104 
STA5H1C5 
STASH! 06 
STASH107 
STASH108 
STASH! 09 
STASH! 10 
STASHlil 
STASH! 12 
STASH113 
STASH! 14 
STASH115 
5W  it  16 
STASH! 17 
STASH! 18 
ST  ASH! 19 
ST ASH 120 
STASH! 21 
STASH! 22 
STASH! 23 
STASH! ?4 


29 


nonn«nrm<?  cm  on  non  non  orct  non  non 
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103 

SO  TO  9999 
faR! T£ I6» l 39) 

104 

GO  TO  9999 

NRI TE ?6» 1 40 ) 

105 

GO  TO  9999 
WRITEI6.I50) 

106 

GO  TO  9999 
kRI TE (6*160} 

107 

GO  TO  9999 
NRITE (6*170) 

108 

GO  TO  9999 
WRITE16,180> 

109 

GO  TO  9999 
WRITEI6*190) 

HO 

GO  TO  9999 
feRITEIGtUOQ) 

ill 

GO  TO  9999 
t)RITE(6»llI0> 

112 

GO  TO  9995 
kRITE(6*U20) 

113 

GO  TO  9999 
URITEI6* 1130) 

50  TO  9999 
5500  5«RlTEt6t5550} 

PHYSICAL  CONSTANTS 

kRITE(6„20f 

«RiTE{6,22JLfK,RHG,CP,0ELTAXIBElTAT,TIHEF 
BOUNDARY  CONDITIONS 
KRI7EI6»30J 

hRSTEI6,33mi,»U2*K2i,X22 
MULTIPLYING  FACTORS 
kRlTEI6,3I/F0,FUFX 
CALCULATION  PARAMETERS 
i»RITE{6r393SNT£RNS»NSTEPX,NSTEPT 
EIGENVALUE  SOLUTION  PARAMETERS 
wRHE  16^61  LI^IT*  ITERHX 


TEST  FOR  TABULAR  CATA 
CO  59  1*1*5 

I FI  NT AS { £  1 *NEeOJ  CAU  VASIN? I) 

CONTINUE 

IF  TABLES  ARE  USeo*TH£¥  ARE  ASSIGNED  AS  FQLLQyS 

sod  functicn  p^iem 

SO.’  FUNCTICN  PHIL?!? 

LO.V  FUNCTICN  PH I PRO  {LANCIA} 

AO. A  FUNCTICN  PH1PF.U  LAH&OA) 

NO. 5  FUNCUCN  FQFXIX} 


5T4SH125 
STASH126 
STASH127 
STASH128 
STASH) 29 
STASHI TO 
STASH1 51 
STASH! 32 
STASHI 33 
STASH! 54 
STASHI 55 
STASHI 56 
STASHI 5? 
STASH! 38 
STASHI 39 
STASHI40 
STASHI4I 
STASHI 42 
STASH! 43 
STASH! 44 
STASH145 
STASHI46 
STASHI4T 
STASH148 
STASHI 49 
STASHI 50 
STASHI 5 1 
STASHI 52 
STASHI 33 
STASHI54 
STASHI 55 
STASHI 56 
STASHI 57 
STASHI 58 
STASH159 
STASH160 
S?^SBl6! 
STASHI62 
STASH163 
STASH164 
STASHI6S 
S7ASH166 
STASHI6? 
STASHI68 
STASHI69 
STASHI 70 

rvs.'tJit  -*  i 
-3  t  ISTv;  f  * 

STASHI 72 
STASHI 73 
S  f ASHI 7< 
STASH! 75 
STASHI 76 
STASH! 77 
STASH178 
STAS  3179 
STASHI 80 
STA5HI81 
STASH182 
STASH! 83 
STASH! 84 
STASH185 
STASHI  8-6 


non  or>  non  non  noon  non 
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STASH! 87 

ERROR  CHECKS  OA  INPUT  DATA  STASH! SB 

ST ASH! 89 
STA5M190 

AUN3£R  OF  INTEGRA! ION  STFPS  MUST  8E  EVEN  OR  ZERO,  ST ASH! 91 

IF  ZERO  THE  PROGRAM  SETS  THE  INTEGRAL  EQUAL  TO  ZERO.  STA5H192 

STASH193 

IFLNSTEPX.EO.Q)  CO  TO  248  ST ASH194 

lFtMGD{NSTEPX«2).EC.O)  GO  TO  248  STASH195 

WRITE  16*5553)  STASH196 

!lRROR=IERROR*1  STASH! 97 


248  IFINSTEPT.EQ.OJ  6C  TO  249 

IF|M00INSTEP?,2).EQ.0)  GO  TO  249 
WRITE{6,5553} 

TERROR* IERRQR+1 


STASH193 

5TASH199 

STASH200 

S7A5H201 


!F( (DELTAX.GT.O. ) 

«  .AND. 

«  lOELTAT.GT.O.I) 

*  GO  TO  250 
WRITE 16*555!) 

! ERROR* TERROR+1 
IFITtMEF.GT.O.)  GO  TO  251 
WRITE (6*5552) 

I ERRORS IERRQR+1 


K12  CANNOT  BE  ZERO  DUE  TO  A  RESTRICTION  IN  THE  FORMULATION, 

IF  U12.NE.0.)  GO  YG  252 
5iRITE<6f5555) 

1 ERROR* i ERROR* l 


CHECK  f=OR  UNDEFINED  SCUNDARV  CONDITIONS. 

252  IFJ ( { K1 1 . N6.0. ) .OR. (K12.NE.0. ) } 

*  .AND. 

»HK21.NE«0.).0R.  (K22.NE.0.)) ) 

»  CO  TO  253 

WRITE (6»44l 
IERRORaIERROR-H 

>53  IF{ IERROR.GT.Qi  STOP 

bRITE(6,333l 

PI*3..14i592& 

ALPHA=K/{RHO«CPi 
WRITEI6*3333)  ALPHA 
CETK*Kn»K22-K12»X21 
WRITE (6*4004) 

WRITEi&,4444)  OETK 

DETERMINE  CONSTANTS  FOR  USE  IN  DO  LOOPS 

CENOH=K12«K22*L-D£TK 

C0HG1=KU*X12«L 

C0N02»K11*«2 

C0N03'=f  {K1?«L )  &«2 

C0N04=Fl 

C0N05~K11«K21 

C0N05=!K12*K22«{  !-**2J  J 

C0N07=»0ETK*L 

CQNQ8=F0 

CO.N09sKI26FL 


STASH202 
STASH203 
STASH204 
STASH20S 
STASH206 
STASH207 
STASH208 
STASH209 
STASH2IO 
STASH211 
STASH212 
STASH213 
STASH2I4 
STASH215 
STASH216 
STASH21? 
STASH2I8 
STASH219 
STASH220 
STASH22I 
STASH222 
5TASH223 
STASH2? 4 
STASH2  5 
STASH226 
STASH227 
STASH228 
STASH229 
STASH230 
STASH23! 
S7ASH232 
STASH233 
STASH234 
STASH235 
STASH236 
:*TASH237 
STASH23S 
STASH239 
STASH24G 
STASH241 
STASH242 
5TASK243 
STASH244 
STASH245 
STASH246 
STASH247 
ST6SK243 


31 
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f.QN10=K22*fO 

STASH249 

CON 11 »F0/KI2 

ST4SH2S0 

CCN12=KI2«DFN0F 

STASH251 

C0N086«K12**2 

STASH2S2 

K ! IK? 1-K1 i *»K2l 

STASH253 

K12X22aKl2*K22 

STASN254 

c 

5TASH255 

#* 

l< 

DETERMINE  EIGENVALUES  FOR  SERIES  SOLUTION 

5TASH256 

c 

STASM257 

IFHK.llK?’-rO.O»  )  *  AN0.IKI2K22.EC .  9.  )  JGO  TO  400 

STASH258 

IFIOET.i  ?r.d#300»100 

STASH259 

100 

CALL  SGI'  tHDETK»EIGENfNTERHS) 

ST4SH260 

GO  TO  900 

STASH2&1 

200 

CALL  SOLVE24  DtTK  »E  IGEN»-NTERHS ) 

STASH262 

GO  TO  900 

SFASH263 

30Q 

CALL  SOLVES !0£TK»EIGEN»NTERHS) 

STASH26A 

GO  TO  900 

STASH265 

4Q0 

CALL  SGLy£4{0ETK»EIGEN»NTERHSl 

STASH266 

C 

STASH267 

C 

SET  UP  INDICES  FOR  00  LOOPS 

STASH260 

c 

STASH269 

900 

NTEHPT*{TIHPF/fiELTAT)*l*5 

STASH270 

M£HPX*{  LABEL  TAX  1+1.5 

STASH2  71 

fcRITE!<*»2flG3} 

ST6SH272 

T*-0ELTAT 

STASH273 

\ 

CO  199  N7=1»NTEKPT 

STA5H274 

\ 

T=T*OELTAT 

STASH275 

\ 

Xa-DELTAX 

STASH2 76 

1 

0099NX«I *  N YEHPX 

STASK277 

6 

IFIT.EQ.0.1  GO  TO  299 

STASH278 

7 

X=X+n£LTAX 

STASH279 

i 

SERIES-O. 

STASH260 

_ 

IF! JPRINT ( D.NE.O)  WRITE! 6»6l00)  T.X 

STASH281 

« 

C 

STASH282 

C 

SET  UP  LOOP  -0  GENERATE  SERIES  SUMMATION  FOR  TRANSIENT  SOLUTION 

STASH283 

' 

C 

STASHZ84 

c 

C0999  l«=l*NTERHS 

STASH285 

f 

ZN=EIG£N! 11 

STASH286 

£ 

9 

C£TAN*ZN/L 

STASH287 

•t 

ZW2*ZN*«2 

STA5H288 

20BN0P*ZN2*0£NCK 

STASH’ 39 

2 

$5NX*$IN18ETAN«XJ 

STASH290 

C8NX=.C0S!BETAN*X) 

STASH291 

SZN=*S  INIZM 

STASH292 

r 

* 

CZN=Ct)S!ZNJ 

3TASH293 

s 

K118N=X1I«BETAN 

STASHES* 

I 

XUBN2=KIIBN*»2 

STASH295 

ZH2SZN*ZN2»SZN 

STASH296 

EXPHUL=EXP{-{ ALPHA# {SET AN) »#2>*r 1 

STASH29? 

SUHCON=BETAN«!KI2«SaNX--KllBN*C8NX)/tZN»!KIIBN2+CON038H-{iai8N2~ 

5T//SH29B 

LC0NO88J »$ZN*CZN— 2.*KIlBM«K12*i SZN++2) ) 

STASH299 

XTERHi®  C {CCNOl*ZN2SZN~CONO?*ZN2SZN*CON03*£  ZN*CZN-SZN)  1 / ZOEnOHj* 

STASHiOO 

IC0N04 

STASH301 

XTERH2®! (CONOS«ZN2SZN<CON06*{ SZN—ZN)+CQN07*ZN*-{ i.-CZN) 1/ZOENOMl* 

ETASK302 

1C0N08 

STASH303 

KFPHIL*KTESI«I*PHIIFQ.} 

SI ASH304 

XFPHIOSBKT6RF2«PHiO!Ovl 

S7ASH305 

XFPHl*EXPFUl » { KFPFIO+KFPHIL) 

STASH306 

1FINSTEPT.EQ.0)  GOTO  399 

STASH3Q7 

TI R|NT»FIAT  (Q.  *T  *NSTEPT  1 7FUHC) 

STASH368 

GO  TO  599 

STASH309 

399 

TIHINT«0. 

5TASH310 

32 
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PI/3, 1415926/ 

^-""slOPE  -0EWtK12K22«L) 

IFCJPRtNT{3).NE.O)  WRtTE(6,6300) 

C0999  IsirNTERPS 
NN*I 

IFCKUK21.EQ.0.)  GO  TO  801 
GO  TO  799 

801  IF{ABS{5LCPE)«LT »1»0)  GO  TO  800 

799  ZU®{2.«HN-i.9*<PI/2.) 

GO  TO  850 

800  Z0MZ.*NN*l.)MPI/2.) 

850  eOUNO»OiO 

CO  99  JslrlTERP* 

U=DeTK»t*Z0/(KUK2l#ZQ9«2+Kl2K22»L**2) 

IF! ABSCSL0PE3  tLT.I.O)  GO  TO  860 
Zl*(NN-I. I«PI+ATANCU) 

GO  TO  870 

860  Z1*NN<*PI+AT£NCUJ 
070  IF { A8StZl-20 }  •  LT .LIMIT)  GO  TO  9 
IFIBOUNO.EQ.cz l-ZO)  GO  TO  90 
BaUNO=Zl“ZO  : 

Z0*ZL 

!FIJPRINT(3).NE.O)  HRTTE<6,«310)  l,J#U#Zl 
99  CONTINUE 

WRITS  56, 11 1  I.ITERMX 
STOP 

9  EiG£Nm=ZI 
GO  TO  999 

90  £'IGEN(I)=CZUZe>72.0 
999  CONTINUE 
RETURN 

ERROR  MESSAGES 


STASH3T3 

STASH374 

STASH375 

STASH376 

STASH377 

STASH378 

STASH379 

STASH3H0 

STASM381 

STASH382 

STASH383 

STASH384 

STASH3B5 

STASH386 

STA5H387 

STASH388 

STASH3S9 

STASH390 

STASH391 

STASH392 

STASH393 

STASK394 

STASH395 

STASH396 

STASH39? 

STASH39C 

STASH399 

STASH4C0 

STASH401 

STASH402 

STASH403 

STASH404 

STASH405 


STASH406 

FORMAT!  lOX<ji2HROGT  NUMBER  ,T3,22HBID  NOT  CONVERGE  AFTER, 15 ,llH  ITESTASH40? 


IRATIONS) 

END  '  - 

i  ■  ■  '  -  .  ' 

IBFTC  S0LV2  M94/2,XR7 

SUBROUTINE  SOLVE  2 1  CEtKp EIGEN f NTERSSS ) 

SOLVES  ?ANUJ=s”DH;*Z/{Xli»K2H*Z*<»2*K22'*Ki2»L*«2r 


SUSH4D8 

5TASH409 

5TASH41Q 

STASH4U 

STASH412 

STASH413 

STA5H414 

STASN415 

STA5H416 

STASH417 

STASH4I8 

STASH419 

STASH420 

STASH421 


Cl HENS ION  EIGEN!  WERHS)  _  STASH* 

£ I KENS ION  JPRINTC3)  STASH417 

OIKENSICN  NTABI51  STASH4I8 

REAL  kUK2l,KI2K22,L,tI«IT,NN  STASH419 

SEAL  LIHIT,L,K,LAH8bA,NUHXSS,NUMS3,KFPHI,KFPHI0,KFPHILfKU,K12  STASH420 
ItX21,K22,N,KTERKl4.KTEfiM2  \  -  STASH421 

COMMON  T,X,LAH80A,KTERMI»KTERM2,ALPHA*BETAi<rfEXPHUL,C8NX,SBNX»DELTASTA5H422 
lX»U€LTAf,Kn,K12,K2I,K22tr-X,NTAff  STASH423 

C0HK0N/R0G?S/NTERHS,L,IT£RHX,XLlK2l,XLZK22,LIfiIT  5TASH424 

CQMHON/PRiNT/JPRiM  -  STASH425 

6300  FORMAT’ IHU46X»24HS0LUT  ION  FOR  ErGENVALUES//26X,7KSOOT  N0fr5X,9HITESTASM426 

•sRATIQN*  10X,  IHU»24X»  IHZ//T  SYASH427 

6310  FQRHATt2SX*I3» 10X,  I3t6X»£15.8, IOZvEIS.S)  STASM428 

CAtA  PI/3.1415926/  STASH429 

JF! JPRIHT!3).NE.O)  WRITE!6»6300T  STASH430 

C0999  I*l*NTEfiPS  STASH431 

J,N«I  STASH432 

ZQ*{NN*l.I*PI_  STASH433 

850  8PUND«0.G  '  STASH434 
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wi  o  nno 


A.FFDL-TR-^-109 


SOI 

799 

@00 

350 


660 

810 


99 


90 

999 


il 


Pl/3. 1415926/ 

SLOPE  *DETK/|Ki2K22*L) 

*  f { JPR I NT  C  3 } *  N£ . 6 }  WRITE {6, 6100) 

C0999  I3i ,NTERP5 
M*l 

IrIKIlSs2t«EQ»0»)  GO  TO  801 
GO  ?G  799 

£HABS{5LCPE).LT.i«0)  GO  TO  800 
20*C2.aNN-1.3«fPI/2.) 

GO  TO  850 

20*C2.*hN*l.)*tPI/2.) 

eOUN0*Q.O 

CO  99  *«I *  I TERPX 

t-=DETK»L*Z0/(KllK2ie20«#2*Kt2K22«L##2) 

IFf A3SCSLCPE1  cLT .1.0)  GO  TO  860 
2l*(NN-l.  I*Pl+ATAN<U) 

GO  TO  870 
21*NN»PI+ATAN{Ul 
IFIA8ST21-20T-l.T-Lf.PTT?  GO  TO  9 
IF I  BOUND. EQ.CZ 1-20 J)  GO  TO  90 
SOUNO*21— 20 
20*21 

1FC JPRINTC3 J.T4E.0)  WRtTEt6,6310)  !,J*U»Z1 
CONTINUE 

NRITCf&tilT  I » ITEGMX 
STOP 

E3GENCX)*21 
GO  TO  999 

IIGENCI)=C2i+ZC)/2.0 

CONTINUE 

RETURN 


ERRCR  MESSAGES 


FORMAT! 10X, 12HRGCT  NUMBER 
IRATIONS) 

END 

O  '  V- 

I8FTC  50LV2  R94/2»XR7  * 

SUBROUTINE  SOLVEZf  CErK.EXGS^NTERSSI 


S7ASU3T3 

STASH374 

STASHT75 

STASH3T6 

STASH377 

STA3H378 

STASH379 

STASH3H0 

STASH38I 

STASH3S2 

STASH383 

STASH384 

STASH385 

STASH3»e 

STASH387 

STASH388 

STASH369 

STASH390 

5TASH391 

STASH39? 

STASH393 

SUSH394 

STASH395 

STASH396 

SX*>H397 

STAS*ri98 

STASH399 

^TASy^vl 

-?'TAS?'i4t5c' 

St$i«402 

5 

szAsmm 


f&samzsz- :s£&i#r&w*nH  i testa zv&tji 

?  5  -  r  ?• .  .  y  sr^SM4os 

■i  yyyy;.;  '.  .  -y  :  stash409 


SOLVES  TANC2J«“09t»Z/(X5t*^t»7«92+A22'»KX29Lf*2} 


STASH4I0 
STASH4I1 
STASH4I2 
5TA5H413 
STASH4 14 
STASN415 
STA3H416 
STASH417 
STASS9I3 
STA5H429 
STSSH420 

stasmazi 


DIMENSION  EIGEMNTEP^tV  --  ' 

DIMENSION  JPRIfcrm  : 

DIMENSION  N TAN {5}  .  :  r  r 

REAL  KUK2I.K12K2 ZtA. 

REAL  t IHI T *1  * K f L 4?§S&sr&3&$£,$!?K53, K F?HI  jKFPHJO, SFP-5 tL S Kil , R l£ 
l  »K21,5<22t  8,  KTERM 

COMMON  T#X*LAM90A/;«i?8iei^§{:S^b>l3L«^S$;rA»f,5XmJLtCB?4S^B^,05S.TASTASHA22 
IX,DELTAttKll,KI2^K?lr52^?«i.JLf.4&  V’  STA5-H4  23 

COHKON/ROOyS/NT^B^LVir.S^T/HmiirSLZR’LfLlNn  S74SH424 

CQPHGh/PRtN?/ JPft&t  ;  SIA5H4Z5 

6300  FORMATJIHl*46X»2NC'iaLT}TX<r^)0OJc  rifc^n VSfc«E5//56X * T^i^T  ftO*5X»9HI TESTASH426 
*RATI0N*I0X*UcUi5«B-t*£//T  .  SYASH427 

6310  F ORHA H 28 X ?  1 3 tL?.S ;* *  S li-ST  STASH420 

GAIA  PI/3.l4i£9?6i<  ;  STASH429 

IFIJPSUNTt3)*«FB»;n  HRTFf'?6^?aO)  STASK430 

C0999  l*s  fNTiXrS  -  .  _T- '  STASH431 

NM*!  STAS K4 32 

20*=CNN*1.  I«PI  STASH433 

8S0  BOO NO =0-0  "  STASH4  >4 
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APFBL-TU-6S-I03 


99 


10 

999 


C 

C 

C 


CO  99  J*l*IT£RMX 

U*OETK#L*ZQ/IK11R21*ZO**24JU2K22*L**2) 
Zl«lNN*PM+ATAMU)  . 
IFIA8Sm-X0J.LT.LlHm  GO  TO  9 
tFteOUNB*£Q.m-Z01:)  GO  TO  90 
B0UND*2I-Z0 
ZO*Zi 

Iff JPRINTI3J »N£.OJ  WRITEI6.6310J- I,J,U.ZI 
CONTINUE 

feftIT£t&,U)  l» JT£RMX 
STOP 

EXGENf I)*Zl  : 

GO  TO  999 

EIGEN M)« (Zl+ZQJ/Z.O 

CONTINUE 

RETURN 

ERROR  MESSAGES 


11 


FPRHATMOXf  12HR00T  NUMBER  ,f3*22H0IB  NOT  CONVERGE  AFTTR.IStUH 
IRATIONS)  - 

END 

X.  '  '  .  . 

SIBFTC  SGLV3  f*94/2,XR^ 

SUBROUTINE  SOLVE3(0£TK,EIGEN,NTERMS-T 
C 
C 
C 


SOLVES  TANTZlsO.O 


6300 

6320 


DIMENSION  EIGEN! NTERMS) 

DIMENSION  UPRIhfl3) 

DIMENSION  NTABI5I 
REAL  K 11K2 1»K 12X22* L.LIMITjNN 

REAL  LIMIT*  L*K*  LAM8DATNUHXS5*NUHS3»KFPMl*K>'PHI0»KFf,HILiiKii*KI2 
ltK2C>K22»N»XTERHl»XTESM2 
COMMON  T9X» LAMBDA* KTERRt*KTERH2» ALPHA* B£TAN»EXPHUL»CBNXr3SNX»0ELTASTASB469 
1 X  »  OELTAT  tK  11 ^  KI2*  K  21 *  K22  *  FX  *  NT  A  B  STASH470 

CQMHDN/RGCTS/NTERKS»L»:tT£RHX,XllK2I,X12K22,LlHIT  STASH471 

COHHON/PRrNT/JPRINT  "c  '  STASH472 

FORMAT!  1H I  *4&3f  »24HSQLUT  T ON  FOR  EIGENVAU>ES//26X,?HR00T  «0,f49X*IHZJ5TASri573 


5USH435 
STA5H436 
STASH437 
5TASI8438 
STASK439 
STASH440 
STASH44I 
5TASJJ442 
STASH443 
STA5H444 
STASH44S 
STA3H446 
STASH447 
5TASH443 
STASH44$ 
STASH450 
STASH451 
STASH452 
STASH453 
I 7ESTASH454 

Stas H4 55 
STASH4S6 
ST4SH457 
SJASH4SB 
ST ASHA 59 
STASH460 
STASH461 
-  5TASH462 
STASH463 
STASK464 
STA5H46<> 
STASH466 
STASB467 
STASH468 


FGRKATI28X, I3,44X,Ei5.S) 

DATA  Fr/3,i4l59?6/ 

Iff JPRINTI3J.NE.6)  HRIT£(6*6300) 

CO  999  r«l,NTERHS 

NN*I 


STASH474 

STASH475 

SJASB4:76 

STASH477 

STASN476 


Zl«NN«PI 

EIGENUJ-Z1  - 

-IFI  JPRINTf  3J.NE.0J  WRltEI S%6320tt1Zl 
999  CONTINUE  :  '  ' 

RETURN  -  ;  ;- 

i  '0  •  v. 

SJ3FTC  S0LV4  7*94/2,  XR7 

SUBROUTINE  S0LVE4f0ETK*EIGEN,NTERKS) 

C 

6  SOLVES  TAMZi=»  INFINITY 
€  - 
DIMENSION  EICENtNIERMS) 

-  DIMENSION  NTABI5J  - 

CIMENSIGN  jpRIfttl3l 
REAL  K11K21,K12K22*L*LIHIT,HN 

REAL  LIMIT a.,X*LAMeDA1.HUMXSS»NUMS5,XFPHIfKF]PHlO»KFPHlL,XiZ,K12 
i  l»R21,K22,N,R*ERHi*KT£RM2 


STASH479 
STASH4B0 
ST4SK481 
STASH4S2 
STASB483 
STASH484 
STASH485 
Sf ASHA 86 
STASK487 
STASH438 
STASK489 
STASH4S0 
STASH49I 
STASH492 
STASH493 
5TASH494 
STASH495 
ST  AS  H4  9  6 


i' 

ir 

t 

| 


I 


i  5 

t  ; 


j 

I 

f 

1 


| 

i 

5 

% 

( 
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AFFDL-TR->(K5-109 


COMMON  T*X,  LAMBDA, KTSRH1  f£TERH2,  ALPHA,  aETAN^ExPXULiS-BN^ijcfcNXt-OELTASTASH^? 
lX»0ELTATiKll»Kl2r»K2l*K22»FXf  NTAB  STASH*  98 

C0BM0N/R0CTS/HTEB«»StUUeRHXJKnK21,K12K22*UH{T  ST  ASH*  99 

CONHON/ERiNT/JPAINT  ST4SH50G 

03C 0  FORMA mHl,46X,24PSGtUTI0N  FOR  EI6£NvAlUES//2sX,7HRaQT  NO, 49X,LH2}STASH501 


6320  FORMAT!  28X,  t3-»44X»El-5*8) 

RATA  Pi/3*t4 15926/ 

TF( iPKlHT {3} ,NE«0>  WRir£(6?630a) 

no  999  i»i,fiT£BKS  - 

KN«f 

n»*<2«*HN-UJ«P.n/2.0 

£IGEMIM*2i 

IF( KJPSlRm J .NE.Ol  HRJTE.T6,«r320Jl*Zl 
999  CONTINUE 
RETURN 
END 

%* 

SI8FIC  SIMPS  P94/2,X*t 

REAL  FUNCTION  FlNTiA,a»NN,F) 

C 
Z 
c 
c 
c 
c 
■t 
c 


SIMPSON  PSTK02  -----  -  : 

slower  acyno 

E*OPPE*  BOUND 

AN-NUK8E8  OF  INTERVALS,  MUST  8fe  EVEN 

F^FUNSHQN  IQ  fig  INTEGRATED,  FUNCTION  HA  HE  HOST  eE  DECLARED 

external  :  -  '  -  .  - 

REAL  LlHI7^t^,LRRffiEiA,^UHXS£,NUHSS^XFPm,fCFPM£0jKFP«lL,Kl2i,A12- 

1 ,  K2.ii%22*h  t  KTERK  l*KT ERH2  _ 

tdhmH  t,X>tAH£0A,KTERKi,KTER«2,ALPHA,86TANe£XPNUL,C8NX,SBNX,DELTASTASH526 


STASH502 
STASH503 
STASH5G4 
STASH505 
STASHS06 
STASH507 
STASM508 
STASH500 
STASH510 
STASH5U 
5 TASKS! 2 
SIASH513 
SfASHSU 
STASH515 
STASK5T6 
STASH51? 
STASH518 
STASHSI9 
STASH520 
STASH521 
STASK522 
5TASH523 
STA5H524 
S7ASHS25 


IX,  DELI  AT,  KO  ,XmB2leK22»FXrNTA8 

m*m:  '  -  -  - 

F=ff»-A!7FA 

SUM^O.^  .  - 

$U%&*QiG  ;  ; 

P»KN-1  .  -  ‘  - 

.  COlOJs 1,8*2  -  - 

.'  .  FJteJ  .  , 

XX-=Fu*H  :  -  -  ' 

10  $y8A*SUHA4FTXX* 

-  .-  PrN.M-2  ^ 

Ca20KK»2*P,2 

FK*KX  ;  :  -  :  - 

XXf'F&H  '  - 

-SUK3*SUK8»FIXXJ  -  -  i 

20  CONTINUE  . 

F*NT={ H/3* 1*IF(A }*f 183*  4„«SUHA>2, *SUH3l 
RETURN  '  - 

END 
S« 

II B FTC  FUMCX  '  P9A/2tXR7 

"Function  xfuhclxxj 
c 
c 
c 


THIS  FUNCTION  sets  UP  THE  INTEGRAND  FOR  Jug  /  INTEGRAL 


STASH527 

STA-SK328 

STASH529 

STASK530 

STASH531 

STA5H532 

STASH53 

STA5B534 

STASH535 

STAS^SJi 

STASH537 

STASN538 

5 TASKS 39 

3TASH540 

STASKS41- 

STASK542 

ST ASKS 43 

STASb544 

STA5H545 

5TASH546 

STA5KS4T 

SF4SHS48 

STASH549 

STASK55Q 

STAS H5 51 

STASK552 

STASH553 

5TASH554 


GIHENSJGN  HIAe*5) 

REAL  lIHITjL,X,LAP0DA,h'UBXSS,NUPSS*XFPHI,KFPHtO,NFPHIL,KlIeXI2 
l,X22,K22,N,KTERBl*KTER82 
COMMON  T,X,LAM8DA,iCTERHUKT£RK2, ALPHA, 8£TAN*£XPNUL,C8NX,S8HX,0EtTASTASHS§5 
lX,3ELTAT,KJl,KI2,X2ItK22,FX,NTA0  -  STASH556 

8NXX*8£TAJ,*XX  :  STASK557 

XFUffC=£XPrUi.»*  KI2*SlNI8NXX3~XIl»6£TA?4*COS*BNXX-J3*FX*FOFXfSXI  7TASH55S 


AFFBL-TXt^-iOS 


hETUSN 

EHD 

TIDFTC  FUHCT  *<H/2,XR7 

rUNCT SCfi  TFUNCftAI'BBAJ 
C 
C 


THIS  FmCtlQH  SFTS  UP  THE  1HTEG8AS0  FOR  THE  LAH50A  INTEGRAL 

REAL  tl«IT»L.«»LAPSDA»?iU}«5S^55UHSS,XfPHl,!a;P{;lO,KPPHIi5ia«,K£? 
l,J%2L*K22oH*JCi'£8HI*KT£ftK2 
COHW3M  r*X,LiKSUA.XT£RKi,KTESH2,A£?HA»8£TAHIE}'PPrULfCE?iXfS3HX,0ElTASrASK569 


ST4SH559 
STfeSHSfcO 
STASH561 
5TASH5&2 
ST4SH563 
$T*5«5fe4 
ST4SH565 
STAS K5 66 
STA5H567 
STASHS&S 


IX»DELTAT*X11 *fa2,K2i*K22*FX,KTAB 
T£SHl=JCT£RHle>PHT?3ifLARaOAJ 
TERK2*xrrAI*2*PPiPRetL*BBD« 

TFUHC=ITERH1^TE?.P2T«£XPT  i ALPHA*! SET AN}**2J  *T  LAKSOft-T  51 
RETURN 
.  END 

S* 

SISFTC  TABLE  P9V2»XS7 
Saa*TOUTIHE  TA81HII1 
C 

c  reads  in  tabular  cat A 


no 

150 

TOO 

50 


fc 

c 

c 


FQRHAfi *5*2F10-O*55H 
-»  T  - 

7=GRHATUHS,61X,9KTA8L£  H0..I2J 
FQRSAT  UHC,  SOX,34HIN0F?£K3SHr 
*8L£  V*RIA8L£1 

FORNATI ASX* £15.8,  H5X, £15*31 
REAL  IKSfAR 
CT«£HSICH-HTABC55 

ClfiESSICX  !?TABtEI51*aHQVAIU5?50t»SEr  Olf5*SOJ 

COSfSJ.WPTSATA/  TNDMR,Q£F_¥AJS 

REaDI5iri)H*4B16Cl> 

kfClT£56»2>  HIieLEM? 

kRIT£l6*3J 

CC  150  J*l*50 

«fEADI 5,lJH*lEDyARII , J J*CE?»ARi 1,U J 
IF3HTABL€m*H>  53*100*110 
feRITEI6o*l  THOVARI 1* J5,CEPVARI I*J1 
<QHELmZ 

fcSIT£?6»5J  STASLECU 

RETURN 

5RIT6f6,6V 

STS? 

ERROR  RESSASE5  - 


STASK5?Q 
STASH571 
STA3K5T2 
STA5H5T3 
ST ASKS 74 
SSASHSTS 
5TASH5T6 
ST ASKS 77 
SKSii>l8 
ST4SKS79 
STASHSCO 
STASK5SI 
STASH5S2 
ST&SB583 
5TASH5S4 

0£?Ef®£HT/5tX»33H  VAX!  AS7ASH535 

ST4SK5^6 


5  F0RS*mHA»41X»9HTA3L€  CCHTA1HS  HGRE  TKJJ4  50  P3IHTS1 

6  FCSHATTlHA»4SXi  27H5RX2R  £?f  TABULAR  INPUT  BA*XS 
EHO 

£»  ~ 

FTSFTC  EHIERP  *9*/2,XRT 

SEAL  FONCTTOS  IflTEFP'TX,!) 


C 

C 


LINEAR  IKTERPOLATICH  FOR  TABULAR  SATA 
CTF-EH5IGS  JTTASISi 

ei££R$J£M  INOVARf 5*507, CSPYA8!  5*501 
COWROS/PTCATAF  IKS¥A8*BEP¥A£ 

SEAL  XKuVAR 
CO  10  J*U50 


STASH537 

STASK5SS 

stashsso 
Sr*S«S9Q 
SSASH5SI 
STAS H5 92 
STA5H593 
STAS H5 94 
STAS5S59S 
STASJ596 
STASS567 
STASH593 
STASHS99 
FLASH* GO 
SjASBsOT 
ST*Sl?S02 
STASHS03 
STASH604 
STASfiSOS 
STASKF06 
5TAS5S07 
STASH603 
5TASH&C9 
ST*SH61S 
STASH5H 
STASH® £2 
STASIS  13 
STA5H&J4 
STASK5X5 
STASN516 
STASH&17 
ST4SKS18 
STASMSi9 
£TA  S3&20 


57 


AWFOI  -TR-6B-109 


1F(TX-INDVAH(L*J)}  20,30, tO 

STASH621 

la 

CONTINUE 

STASH62? 

UK)TE16,99)  i 

STASH623 

CALL  EXIT 

5TASH62A 

20 

INTERP*D£PVAfmt  J-1)*{CEPVAR(  1 ,  S-DEPVARS 1 ,  J-15  )  I TX-INDVAR U » J~ 

1 )STASHfc2S 

*if{  'NDVARt  l»J5-IN0¥AfU  1,  J-lll 

STASH&26 

GO  TO  100 

STASH627 

10 

IN?£R.P*DEPVAR«  I,  Jl 

STASH62B 

SOG 

RETURN 

STA3H629 

C 

STASH630 

w 

EKRCR  MESSAGES 

STASHA11 

c 

STASH6  >.? 

99 

FORMA T { 10X, 16HARGUMENT  EXCEEDS  EXTENT  OF  TABLE  NO. ,12) 

ST  ASHA  11 

END 

STASM6 14 

E# 

STASH61S 

£SBf 

TC  PH  10 T  094/3, XH? 

STASH616 

FUNCTICN  PHIOSTT) 

STASH617 

C 

STASH618 

j- 

THIS  FUNCTION  CALCULATES  THE  INSTANTANEOUS  VAt-UE  OF  THE 

STASH&19 

C 

TINE-VARYING  BOUNDARY  CONDITION  AT  X*0. 

STA5H&40 

c 

fHfc  FUNCTION  Pf-IOfT)  NAY  BE  tOACFD  INTO  THE  PROGRAM 

STASH6AI 

c 

AS  AN  ANALYTICAL  EXPRESSION  OR  AS  POINT  DATA  IN 

STASH642 

c 

TABULAR  FCRM 

STASH641 

c 

STASH6AA 

REAL  1NTERP 

STASH6A5 

REAL  LIMIT,L,K, LAMBDA, NUHXSS,NUNS3,KFPHI,KFPHI0tKFPHIL,K 11 »K12 

STASHf.46 

l,K2ic.R22»N,XTefl?4l,KTE«N? 

STASH64? 

COMMON  TtX,lAHeoAfKTE«MI,KU'RM2, ALPHA, BET AM, EXPHUL,C8N\,S3NX, DELTAS! ASH648 

lX,0ELTAT,Kn,*U2,K2I,K22,FX?NrAS 

STASHBA9 

DIMENSION  NTABS5) 

STASH6S0 

IFCNf  AB{ I I.NE.O)  GO  TO  100 

STASH&9I 

c 

STASHSS2 

c 

PH 10= ANY  FUNCTION  OF  TIME 

STASH651 

c 

STASH654 

PH  IDs.  1.0 

STASH6S5 

RETURN 

STASH656 

100 

PhIO=lNTERP(TTtl) 

STASH&57 

RETURN 

STASM6S8 

END 

STASH659 

$» 

STASH66C 

$I8FTC  PHUT  P94/2.XR7 

STAS H&6L 

FUNCTION  PHIL  ITT) 

"T4SH 662 

C 

STASH661 

c 

THIS  FUNCTION  CALCULATES  THE  INSTANTANEOUS  VALUE  OF  THE 

STASH&6A 

r. 

TIME-VARYING  BCUNCARV  CONDITION  U\  X«L. 

STASH665 

t 

THE  FUNCTION  PHI  LIT  5  MAY  BE  LOADED  INTO  THE  PROGRAM 

STASH&U6 

c 

AS  AN  ANAtYT ICAi  EXPRESSION  OR  AS  POINT  DATA  IN 

STASH6S7 

c 

TAJULAR  ECR*. 

STASH* 6 8 

c 

STASH6&9 

REAL  INTER? 

STASH670 

REAL  LIMIT, L,K,lSMEDAfNUMASS,NUMSS,XFPHI,KEPHI0,KFPHIL,KlI,X12 

STASH&  ?1 

i,K21,k22,K,XTERP  UKTERN? 

STASH-S72 

COMMON  T , XstSHBC  A* XT  JVM I »K7ERH2, ALPHA, 8ETAN,CXPMUL ,C0NX,SBNX,QfiLTASTA$H6  ?1 

!S,DEITAT»S1I,K12,K21,<C2?,FT,NT!:B 

ST ASHA 74 

DIMENSION  NTABI5 ) 

STASHS75 

i Fl NT AB 1 2 ) .Nt *G|  GO  TO  iUO 

STASH*  Jf> 

STASH* 77 

c 

PHIL-ANY  FUNCTION  OF  '  •' 

ST &$K6 1H 

c 

STAS KG  T9 

PHiI.=  I . 

STASH680 

RETURN 

STASM631 

100  PHUa  lNT«-«PlTT,2,i  STASH66? 
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RE TURN 
FNP 

sibftc  qerivo  fva/^.xr? 

FUNCTfriNFHIPROl  LAMBDA  5 

c 

C  THIS  FUNCTION  CALCULATES  THE  INSTANTANEOUS  VALUE  OF  THE 

C  CERIVATIVE  OF  THE  TINE— VARY  I  Ml  BOUNDARY  CONDITION  AT  X*0. 

C  THE  FUNCTION  HAY  8£  LOADED  ANALYTIC ALLY  OR  AS  POINT 
C  OATA  IN  TABULAR  FCRH. 

c 

REAL  INTERP 

REAL  LIHI  I,tfK,LAFBDA»NUHXSStNUHS^,KFPHS*KFPHIOIKFPHH  ,K11,KL2 
l.m,K22,NsKTFRHl,KTEftH2 


ST  ASHA  S'} 
ST  ASHA HA 
STASH&HS 
STASH686 
STASH68? 
STASKAB8 
STASH689 
STASH690 
STASH69I 
SYASH&92 
5TASH693 
STASH694 
STASH&95 
5TASH696 

CGHHOfC  T, X»LAH0OA, KTE8Kl,KTE«H2f ALPHA, B£TAN,EXPHUL*C 8NX, SUNS, 0FLTAS1ASH69? 
lX»0cLTAT»KIl,K12,K?I,K22*FXfNTAB 
01 HENS I ON  NTABISJ 
IFINTAB531  ■»NE.Q)  60  TO  iOO 


C 

c 

c 


100 


PHIPRG=ANY  FUNCTION  QF  T i  “T£ 

PHIPR0«O. 

RETURN 

PHI PRO- INTf RP{ LAHBOA»  3} 

RETURN 

END 


*  =■ 

6I8FIC  DERIVL  H9A/2»X«7 

functicnphiprulambdai 
c 
c 
c 
c 
c 
c 


THIS  FUNCTION  CALCULATES  THE  INSTANTANEOUS  VALUE  OF  THE 
DERIVATIVE  OF  THE  TIHE-VARYING  BOUNDARY  CONDITIONS  AT  X^L* 

THE  FUNCTION  HAY  BE  LOADED  ANALYTICALLY  OR  AS  POINT  DATA 
IN  TABULAR  FORK 

REAL  IN7FRP 

PEAL  I.IHI  T  t  >.»K,LAf'BOA*NUHXSS*NUMSS»KFPHIvKFPHIO»(<FPKIL»XH  *Ri2 
itK2M;,;2?N»KTERNlfKTEftH2 
COWGN  T, X*LA.;30A-KTERHl?KT£RH2f  AI  PHA»BETAN,EXPHUL,CBNX,SSNX,OEci  ASTASH721 


STASH698 
S1ASH699 
STASK700 
STASH.  J1 
STASH? 02 
STASH703 
STASH? 04 
STASH705 
STASH706 
STASH707 
SI4SHY08 
STASH709 
STASH710 
STASH7H 
STASH7L2 
STASH713 
STASH71A 
STASH715 
STASH? 16 
STASH71? 
STASH7I8 
ST  A$!57  19 
STASH720 


iKtOELTAT  I  ,K  •  !  vK??*FX»NTAB 


TC 


C 

c 

c 


A) 


♦I8FTC 

C 

c 
c 

c 
c 


STASH? 22 
STASH723 
STASH724 
STASH725 
STASH726 
STASH727 
STASH728 
STASH729 
STASH? 30 
STASH731 
STASH732 
STASH733 
STASH734 
5TASH735 
STASH73& 
STASH? 37 
STASH738 
STASH739 
STASH7A0 
S  ASH?^I 
Si 4SH742 
STASH7A3 

T»X5LAMsOA,KTfH«l  r-u.  ,  ALFHAfc8ETAN» EXf’HUl  rCSMX-SfiNX, DELTAS!  ASH7&A 


If  INfABKKNL.OI 

PHJPRL=ANt  FIJ-'CT,  h 

FHlfRLsG. 

RETURN 

LOO  PHIPRi* INTERPILANBDA 
RETURN 
END 
l* 

FX  P94/2*XS? 

UNCTICNFCfXi...  5 


THIS  FUNCTION  CONFUSES  THE  |*Jl»fAl  CONDITIONS  OF  THE  inO, 

TKLSE  INlTiAl  CONDITIONS  S*Y  BE  tOAOED  INTO  THE  PROSRAn 
ANALYTICALLY  oa  AS  PQINf  OAf*  IN  TABULAR  TORN. 

RFAl.  1NTERP 

RfAL  LIH!  (,L»KeLA?'BDA»NUHEASlNuKrv5»KFP»{*»KFPMlO*KrPNtl#xn  ,£  i2 
1 7K?.i,r. 

COMHOS 


«  w  o 


1.X»mTA¥«XllTKi2*K2if»22tf 
CJHSNSIGN  H?Ag{5J 
IHN?A8(5}*NF,0}  GO  TO  ISO 

FOPS*  mi  FUNCTION  Of  X 

FOfXm, 

5$  f?  £  TUSH 

IQS)  POF'<*tNTERP(SX?5l 

am**# 

CHS 


3rASH?45 
STA5H746 
S?A!>H74? 
SU5H748 
STASH?*'} 
STASH750 
STA5H751 
$T4SH7« 
STASH?  S3 
STASH7S* 
STASH7SS 
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APPENDIX  Tl 

EIGEN V  ALIJS  SUBROUTINES 


The  solution  of  equation  S3  depends  on  values  of  f3n  and  which  are  derived  from  the  positive 
eigenvalues  of  F  nation  40, 


DU, 


t«  n 


iR  '  Ks,K„2Z  +  KS£K)2  L* 


(40) 


-Since  we  are  seeking  positive  valuer  of  the  sign  of  tbs  left-hand  side  of  the  equation  may 


he  associated  with  tbs  parameter ,  D.  Thus  three  formulations  are  possible  corresponding  u> 
Es  being  positive,  negative  or  zero.  A  fourth  possibility  is  that  of  the  denominator  going  to  rc.ro. 
The  last  two  solutions  are  trivial.  If  D  is  zero  we  have 


tan  zn  =  Q 

Tne  solution  to  this  equation  is  merely 

zn  =  n  fr 

If  the  denominator  of  equation  40  goes  to  zero  we  have 


(54) 


The  solution  to  the  equation  is 


ton  zn  -  a 
(2n  ~  I )ir 


2n  = 


i5) 


However,  if  D  has  a  value  other  than  zero  the  equations  are  solved  by  an  iterative  process. 
In  these  cases  the  eigenvalue  subroutine  has  been  programmed  to  ignore  the  root  a 1  (0.0)  since 
this  produces  a  trivial  solution.  The  procedure  then  is  'outlined  below  for  a  positive  value  of  D. 


We  know  that  the  solution  lies  between 

r  (2n  "  !  )  T 


{ n  -  I )  tr, 


n  =  i ,  2 , 3, 


For  a  first  approximation  to  the  root*  zq  q,  we  shall  choose 


We  then  write  two  equations 

u 


{gn  ~  I )  7f 
2 

D  L  Z, 


.m-i 


n,  m 


^21^11  ^,m-l  '  ^12 


a,  m 


=  (n-i)  u  -t  tan 

,  ■  th  . 

refer  to  the  m  it 

then  determined  to  any  desired  LIMIT  of  accuracy  L*y  writing 


^  uo, 


( 56) 

(5') 

(58) 


Where  the  subscripts  n  and  m  refer  to  the  iteration  toward  the  n^1  root.  The  root  is 


** 
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Provision  is  made  in  the 
The  oasic  difference  in  the  5 
zn,o’  For  negative  values  of 


•iud  the  solution  proceeds  as 


^  <  LIMIT 


( 59} 


solution  for  a  negative  ^ue^of^D^r  Ojea  ^ n 
D  Equation  56  becomes  R  approximation. 


~  Cn  +  I )  w 

before  with  Bquatioi  58  becoming 


(60) 


=  n  ^  *»ai» -*{,,  , 


(61) 
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TABLE  IT 

COMPARISON  OF  DATA  FOR  CASE  1 


T(0  ,t)  =  t ;  T(L,t)«0;  T(x,0H> 


TEH?  3 

TEH?  @ 

t=»50 

TEMP  3 

t-500 

TEMP 

t=1000 

L 

£=0) 

ST  iSH 

LTA 

STASK 

LTA 

STASH 

LTA 

o.,o 

0, 

00 

50,00 

50.00 

500.00 

500.00 

1000.00 

1000,00 

0,5 

0. 

00 

37.39 

37.38 

455.95 

455.93 

930.71 

930.09 

1,0 

27.45 

27.46 

414.85 

414.83 

864.38 

864.35 

1*5 

19.79 

19.80 

376.56 

376,53 

300,87 

$00.83 

1,0 

13.99 

14.00 

340.92 

340,87 

740.02 

739.97 

2,5 

9.76 

9.70 

307.74 

307.69 

681 .66 

681.60 

3.0 

6.58 

6,59 

276.88 

276.83 

625.65 

625.59 

3.5 

4.37 

4.38 

248.13 

248.12 

571.82 

571.76 

4.0 

2  .84 

2  85 

221.48 

22; ,42 

20.03 

519.96 

4.5 

1.80 

1.81 

196.61 

196.56 

470.11 

470.0a 

5.0 

1.12 

1.13 

373.43 

173.37 

42.1.90 

421.84 

5.5 

0.68 

0.68 

Jt  51.77 

151.72 

375.26 

375.20 

6.0 

0.40 

O.AJ 

131.43 

131.43 

330.03 

329.97 

6.5 

0.23 

P.,24 

132.44 

112.36 

286.04 

235.93 

7.0 

0.13 

0.33 

94.38 

94 . 34 

243.15 

243.10 

7.5 

e.ov 

0.07 

'7.27 

77.24 

201.19 

201 . 15 

8.0 

0.04 

1 0-91 

60.89 

160.02 

159.98 

8,5 

£.02 

0.02 

45.16 

45.14 

119.67 

119.44 

9.0 

0.01 

0.01 

29.86 

29.84 

79. 3S 

79.37 

9.5 

0.00 

0.00 

14.85 

14.84 

3^,61 

39.60 

10.0 

0. 

00 

0.00 

0.00 

0.00 

L_ 

0.00 

0.00 

0.00 
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TABLE  III 


COMPARISON  OF  BATA  FOR  CASE  2 


T(Q,t)*=lCO+t;  T<L,t)«102j  T(x,C)»102 

: 

TEHP  Q 

TEMP  ( 

1  t*»50 

TEMP  @  t^SOO 

TEMP  @ 

t=10QO 

L 

t=0 

STASH 

LTA 

STASH 

LTA 

STASH  j 

LTA 

0.0 

100.00 

150. GO 

150.00 

600.00 

600.00 

ilGO.OO  1 

1100.00 

0.5 

100 

.00 

137.33 

137.38 

555.95 

555.51 

1030.71 

1030.69 

1.0 

127.45 

127.46 

514.85  ' 

514.83 

964.38 

984.35 

i.5 

119.79 

119.80 

476.56 

476.53 

900.87 

900.83  ' 

113.99 

114.00 

440,92 

440.37 

840.02 

839,97 

2.5 

109.70 

109.70 

407.74 

407.69 

781.65 

731 .60 

106.58 

106.59 

376.88 

3/6.83 

725.65 

725.52 

3.5 

104.37 

104.38 

348.18 

348.12 

671.82 

671.76 

4.G 

102.34 

102.85 

321.48 

321 ,42 

620,03 

619.96 

4.5 

101 .80 

101.81 

296.61 

296 ,56 

570.11 

570.04 

5.0 

101.11 

101.13 

273.43 

273.37 

521.40 

521 „ 84 

5.5 

100.68 

100.68 

251.77 

251.72 

475.76 

475. 2U 

6.0 

100.40 

100.40 

231.48 

232.43 

430.03 

429.97 

6.5 

100.23 

100.24 

212.40 

212.36 

386.04 

385.9? 

7.C 

100.13 

100.13 

194.38 

194.34 

343.15 

343.10 

7.5 

100.07 

100.07 

177.27 

177.24 

301.19 

301  a5 

8.0 

100-.04 

-100.04 

160.91 

160.88 

260.02 

259.98 

8.5 

100.02 

100.02 

145.  .6 

145.14 

219.47 

219.44 

y.o 

100.01 

100.01 

123.36 

129,84 

179.38 

179.37 

9,5 

t 

100.00 

100.00 

114.85 

114.85 

139.61 

139.60 

ioo .eo 

100.00 

100,00 

100.00 

. 

100 .OC 

10G.00 

100.00 

_ 

Z 

_ 
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TABLE  iv 


COMPARISON  OF  BATA  FOR  CASE  3 


- 

T(0,t)=t;  T(L,t)~102;  X<x,0)« 

' 

TEKP  @ 

TEH?  @ 

£*50 

TEMP  Q 

t=59G 

tek?  q 

t*ioeo 

L 

t«0 

STASH 

LTA 

STASH 

LTA 

STASH 

LTA 

0.0 

0.00 

50.00 

50.00 

500.00 

ir  j 

SCO. 00 

1000.20 

1000.00 

0.5 

5.00 

42.38 

<2.33 

460.95 

460.93 

935.71 

S35.69 

1.0 

10. 30 

37.46 

37.46 

434.35 

424.83 

874.38 

874.35 

1.5 

15.00 

34.73 

34,80 

391.56 

391,53 

815.87 

815,83 

2.0 

20.00 

33.99 

34.00 

360.92 

360.87 

760.02 

759.97 

2.5 

25.00 

34.70 

34.70 

'32.74 

332.69 

706.66 

706.60 

3.0 

30.00 

36.58 

36.59 

306.88 

306.83 

655.65 

655.59 

3.5 

35.00 

39.3? 

39.35 

283.18 

283.16 

606.82 

606.76 

4.0 

40.00 

42.84 

42.85 

261.48 

261.42 

560.03 

559.96 

4.5 

45.00 

46.80 

46.61 

241.61 

241.56 

515.11 

515.04 

5.0 

50.00 

5L.1I 

51,13 

223.43 

223.37 

471.90 

471.84 

5.5 

55.00 

55.68 

55.68 

206.77 

206.72 

430.26 

430,20 

6.0 

60.00 

60.40 

60.41 

191.48 

191.43 

390.03 

359.97 

6.5 

65.00 

65.23 

65.24 

177.40 

177.36 

351.04 

350.99 

7.0 

70.00 

79.13 

70.13 

164.38 

164.34 

323.15 

313,20 

7.3 

75.00 

.  75.0/ 

75.07 

152.27 

152.24 

276.19 

276.15 

S.O 

80.00 

SQ.C's 

80.04 

140.91 

140.89 

249.02 

239.98 

8.5 

85.00  -1 

:85.02 

85.02 

139.16 

130.14 

204.47 

204.44 

9.0  „ 

90>30 

90.01 

90.01 

119.86 

119.84 

169.33 

169.37 

9.5 

95.00 

95.00 

95.00 

109.85 

109.84  i 

134.61 

134.60 

10.0 

100.00 

1C0.00 

100.00 

109.00 

100.00 

ICO .GO 

100,00 

48 
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&- 


CGMPAF. 


1  ij-.-f. 


TEMP  @ 

TEH?  7 

£=50 

TEMp  S 

i  =500 

* 

-190- 

* 

L 

t=U 

STASES 

LTA 

STASH 

_  LT*  ’ 

7/A-,  { 

iiI7.  -  - 

* 

0.0 

0..0 

50.00 

50.00 

500 .00 

5Sa  .0-.- 

~ 

ir&j.m 

i 

i 

0.5 

5.00 

42.3-3 

42.38 

4-50.95 

460.5- 

535.71 

S3*.,  J9 

i 

* 

i.o 

10, OO 

37.46 

37.46 

434 . S5 

424.33 

5.*  - ,  3d  i 

874.'"- 

i 

1.5 

15.00 

34.79 

34.80 

391.56 

391.53 

£15. -3/ 

£15,63 

* 

2.0 

20.00 

33.99 

34.00 

360.92 

360.87 

760 .02 

7  3 1* .  V  7 

i 

2.5 

25.00 

34.70 

34.70 

332.74 

332,65 

706.66 

706. 6C- 

i 

3.0 

30.00 

36.53 

36 . 50 

306.88 

306.33 

655.65 

665.59 

i 

/»  r* 

J 

35.00 

39.37 

39.38 

233.18 

283.1o 

606.82 

506,76 

i 

4.0 

40.00 

42.34 

42.85 

261 ,4S 

261.42 

56Q.03 

559.96 

i 

j 

A. 5 

45.00 

46.80 

4n4bi 

241.61 

241.56 

515.11 

515.04 

5.0 

50.00 

5iai 

51.13 

223.43 

223.37 

471.90 

471. S4 

5.5 

55.00 

55.63 

55.68 

206.77 

206.72 

430.26 

430.20 

6.0 

60.00 

60.40 

60.4} 

191.48 

191.43 

350.03 

389.97 

6.5 

65.00 

.5.23 

65.24 

177.40 

1/7.36 

351.04 

350.99 

7.0 

70  .eo 

70.13 

70.13 

164,38 

164.34 

313.15 

313.10 

i 

7.5 

75.00 

75.07 

75,07 

152.27 

152.24 

278.19 

276.15 

j 

g.O 

SO  .00 

80.04 

80.94 

140.91 

140.89 

240.02 

235.98 

3.5 

S5.0G 

85.02 

85.02 

130.16 

00.14 

204.47 

204.44 

l 

9.0 

90.00 

90.01 

90.01 

U9, 86 

115.84 

189.38 

169.37 

i 

i 

9.5 

>i»00 

95.00 

95.00 

I 09. S3 

109.84 

13*.61 

134.60 

i 

10.0 

103.00 

ICO. 00 

100.00 

100.00 

100.0-0 

100.00 

IOC .09 

i 

_j 

Or* 


V*.  \ 


fc:.  i 
L  c  i'4 


^  ^  d 
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TABLE  V 

CAPARISON  OF  MXA  FOE  CASE  4 


T<0>e)*2000;  q<L,t)-0|  !{*,0)=1G9  | 

TEMP  @ 

TEMP  $ 

£*300 

C-1SG0  ~ 

TEMP  @ 

t-3600 

L 

LTA 

STASH 

LTA 

STASH  | 

LTA 

l"  “1 

' 

100.60 

1000. GO 

iOOG.OO 

1000.00 

1000.00 

, 

1300.00 

1000.09 

.0 

100.00 

S59.95 

899.90 

994-85 

994.40 

999.85 

999.80 

2.C 

100.00 

£02.53 

802.50 

989.83 

989.45 

999.71 

999.65 

3,0 

100.00 

710.26 

709.99 

985. 0b 

984.72 

999.57 

999.53 

4.0 

625.48 

625.30 

980.06 

978.98 

999.44 

999.41 

5.0 

■  ■ 

550.05 

976.72 

976.13 

399.33 

9*9.2? 

6,0 

486.38 : 

486.14 

073.37 

972.66 

999.24 

399,19 

7.0 

435.25 

434.97 

970.67 

970.04 

$99.16 

999.12 

8.0 

160-00 

593 .98 

**0*7  CO 
>  #»r*- 

968.08 

999.10 

999.06 

9.0 

100.00 

375.51 

375.02 

967.49 

967.10 

999.07 

999.91 

10.0 

100. GO 

367.70 

367.56 

L _ _ 

967.03 

. 

966.99 

999.05 

' 

598.99  [ 

TABLE  VI 


COMPARISON  OF  DATA  FOR  CASE  5 


T(0,t}=0;  q{i, 

»*)"  “sec5  t'x»0Ja° 

7 

- n 

TEMP  § 

TEMP  0 

e-300 

T2*P  fe 

t-1800 

TEMP  § 

£-3600 

L 

$Xa9b 

LTA - 

STASH 

LTA 

STASH 

LTA 

0.0 

0.0 

0.0 

O.t 

G.O 

0,0 

0.0 

0,0 

1.0 

0.0 
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Figure  10.  Temperature  Profiles  (Casa  3) 
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Figure  12.  Temperature  Profiles  (Case  5) 


